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ABSTRACT 


This document Is a description of and user’s manual for 
one of a group of FORTRAN programs which/ together/ can 
be used for the analysis and design of wings In steady/ 
subsonic flow according to a kernel function method 
lifting surface theory* This particular program Is the 
one which solves the sets of simultaneous/ linear/ 
algebraic equations arising from the thin wing 
analysis. This program has the capability of striking 
out rows and columns of the aerodynamic Influence 
matrix and rows of the associated boundary condition 
vectors (right hand sides). This capability 
significantly enhances the effectiveness of the kernel 
function method of lifting surface theory because 
studies of the convergence of solutions with the number 
of control points can be done with the calculation of 
only a single Influence matrix. 
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Richard T. Medan 
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and 

Opal J, Lemmer 

Computer Sciences Corporation 


1 INTRODUCTION 


This document Is a description of and user’s manual for 
one of a group of FORTRAN programs which/ together/ can 
be used for the analysis and design of wings In steady/ 
subsonic flow according to a kernel function method 
lifting surface theory. This particular program Is the 
one which solves the sets of simultaneous/ linear/ 
algebraic equations arising from the thin wing 
analysis. This program has the capability of striking 
out rows and columns of the aerodynamic Influence 
matrix and rows of the associated boundary condition 
vectors (right hand sides). This capability 
significantly enhances the effectiveness of the kernel 
function method of lifting surface theory because 
studies of the convergence of solutions with the number 
of control points can be done with the calculation of 
only a single Influence matrix. 


The theory behind the program Is described In ref. 1 
while the other/ associated computer programs are 
described In refs, 2-6. 


Questions concerning either this document or the 
computer program or the associated computer programs 
should be directed to: 

R. T. Medan 
Mall Stop 247-1 
Ames Research Center 
Moffett Field/ 

Ca. 94035 
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2 GENERAL QESQRLRtmM 


This computer program solves sets of simultaneous/ 
linear equations derived from an aerodynamic Influence 
matrix which Is stored on disk along with the boundary 
condition matrix. After the solution Is obtained it is 
stored on a disk file for later use by programs which 
eva.luate loads, pressure distributions, etc. 


The program can be operated either In a batch or 
conversational mode. if done conversationally the user 
Is prompted for Input. By entering various commands 
the user can have the solutions stored or printed, the 
boundary condition matrix printed, spanwlse and 
chordwlse control points printed, etc. Operating 
Instructions are given In section 4 while a complete 
description of the commands Is Included In sections 5 
and 5. 


The program has the capability of deleting rows and 
columns by the use of the DSM, DCM, DSP, DCP, and DPP 
commands. Deleting rows deletes control points while 
deleting columns deletes pressure modes. After a case 
Is completed the user has the option of: (1) going 
through the same file again deleting more rows and/or 
columns, (2) starting a new wing, or (3) halting the 
program. 


If the wing Is symmetric, then the program will employ 
symmetry to reduce the amount of computation. In the 
case of wing symmetry the boundary conditions stored on 
the boundary condition file can be symmetric, 
antisymmetric, or both and the aerodynamic Influence 
matrix can contain symmetric, antisymmetric or both 
types of spanwlse pressure modes. If the boundary 
condition file contains both types of boundary 
conditions and the Influence matrix file both types of 
spanwlse modes, then the user can obtain solutions for 
symmetric cases only, antisymmetric cases only, or 
both. 
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3 METHOD OFSQLUXLQN 


Householder’s method for solving simultaneous linear 
equations Is utilized. This method was chosen because 
it is applicable to both square and rectangular 
matrices and In the case of rectangular matrices It Is 
not necessary to least square the matrices first (ref. 
7). 


There is a storage advantage In that the whole matrix 
need not be In core at one time. The program reads in 
only as many rows as there are chordwise control 
points. The method was taken from ref. 7 and modified 
to achieve the capability of striking out rows and 
col umns . 


This subroutine t r i angul ar I zes the influence matrix by 
means of orthogonal transformatl on matrices which 
preserve the conditioning of the matrix. ThIS/ along 
with the elimination of the need for a least squares 
routine greatly improves the numerical accuracy and 
stability of the solution over that of the Gaussian 
reduction method. 
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4 ua£Rs maTjmcTioNs 


4.1 LfiLILAL S.ETU.P.-AMES V Jjaoy.SMSS 

For either batch or conversational processing the 
following TSS commands must be given. These commands 
are required once and only once for each user ID, The 
first three commands create the Identification number 
file named IDFILE* This file contains four zeroes In 
binary form, 

SHARE MEDAN^FSARTM^INIDFILE 
CDS MEDANJDFILE 
DELETE MEDAN 

SHARE MEDAN^FSARTM^SPROG.Vl 

4.2 ,£QNV.£RSATIQNAl. MQJIE— AMES^ 360/67 TSS 


See section 4.1 for the Initial setup. All Integer 
data should be entered In a 1615 format^ all 
alphanumeric data left justified beginning In column 1. 


After logging on proceed as follows: 


User ! ames usysl I b 

jobl lbs sysul I b 
Jblb medan 

User : ca 1 1 sol vr$ 

Prog: ENTER BATCH 

User: Carriage return 

Prog: ENTER ODISK 

User: For terminal output enter carriage return. 

For output to a disk file enter a positive 
non-zero number less than 10, For the AMES* 
TSS system the output will be found on the 
file named OUTPUT. SOLVR, NX where X Is the 
numerical value of ODISK, The program Issues 
Its own DDEF commands so no control cards are 
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Prog: 
User : 


Prog: 
User : 


neeHeH, The program uses logical unit 4 for 

this output, 

ENTER IDl^ ID2^ ID3 

Enter Identification numbers 

IDl Identification number of the geometry 
file from which the AIM and BC flies 
have been derived. 

ID2 Identification number of AIM file, 

ID5 Identification number of the boundary 
condl t Ion file. 

These Identification numbers are used to 
create DDEF*s which are Issued by the 
subroutines AIMFIL^ BCFIL and BNKFIL. 
This makes It unnecessary for the user 
to Issue DDEF commands on the TSS 
system. These numbers are also checked 
against the Identification numbers on 
the files. This Is unnecessary in the 
TSS version/ but would be helpful as an 
aid In Insuring that the correct files 
are being used If the program were to be 
converted to another computer. 

ENTER NTYPE 

(NTYPE will not be requested for an 

unsymmetric wing) 

Enter : 

-1 to obtain solutions to the antisymmetric 
cases only/ 

0 to obtain solutions to the symmetric 
cases only/ or 

+1 to obtain solutions for both types of 
cases. 
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Prog: 

User: 


At this point the user must begin entering 
commands. After each command he will be 
prompted for supplemental Input or with 
another + sign Indicating that he should 
enter another command. The full list of 
commands Is given In section 6 while a basic 
list for the beginner Is Included In section 
5. After the NEW or KNEW command the next 
Input will be ODISK followed by IDl^ 102/ 
103/ and NTYPE. Use the HALT command to 
terminate execution. 
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4.3 MKE mode- am ES^ 360/6 7 tSS 


See section 4.1 for the Initial setup. All Integer 
data Should be entered In a 1615 format/ all 
alphanumeric data left justified beginning In column 1. 


CARD 1. . . . AMES USYSLI B 
CARD 2....JOBLIBS SYSULIB 
CARD 3. . . . JBLB MEDAN 
CARD 4.... CALL SOLVR$ 

CARD 5....T 

CARD 6. , , , F I 1 e Identification numbers-- 101/ 1D2/ ID3 

IDl Identification number of the geometry 
file from which the AIM and BC files 
have been derived. 

ID2 Identification number for AIM file. 

ID3 Identification number of BC file. 


These identification numbers are used to 
create DDEF*s which are Issued by the 
subroutines AIMFIL/ BCFIL and BNKFIL. 
This makes It unnecessary to Include 
control cards on the TSS system. These 
numbers are also checked against the 
Identification numbers on the files. 
This Is unnecessary In the TSS version/ 
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Card 7. 


NEXT 


but would be helpful as an aid In 
Insuring that the correct files are 
being used If the program were to be 
converted to another computer. 


. ..NTYPE (omit for an unsymmetrlc wing). 

Enter : 

-1 to obtain solutions to the antisymmetric 
cases only. 

0 to obtain solutions to the symmetric 
cases only. 

+1 to obtain solutions for both types of 
cases. 


Follow by cards containing commands 
(described In sections 5 and 6) In columns 
1-3 each followed by a card containing 
command arguments (If the command calls for 
arguments) as specified In section 6. 
Following a NEW or KNEW command Input resumes 
with card 6, The HALT command must be the 
last card of the Input deck. 
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4.4 OTHER COMPUTERS 

Remove all calls to BCFIL^ AIMFIL^ BNKFIL, OBEY, and 
CVRT In the main program and use appropriate tape or 
disk control cards In their place. These, hopeful 1 y, 
are the only changes that need to be made since 
considerable effort was made to code the program in 
standard FORTRAN, Then follow the Instructions, where 
appropriate. In section 4,2 and 4,3, 
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5 BASIC COMMANDS 


These commands will enable a new user to run the 
program Without deleting rows or columns or having any 
options. The solutions will be stored. There Is no 
Input associated with these commands. 


HALT. ..... TermI nates program. This command Is needed 

and Is the last I tern of Input. 


NEW. ...... Beg I n a new wing. 

START Begin processing. 
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6 PULL PEaCRIF-ILOH QL CQMMAfiP-^ 


In this section the commands which control the flow of 
the program are described* In all cases the first 
three letters of a command are sufficient Input. All 
numeric fields are to be entered In a 1615 format. All 
alphanumeric data Is to be left justified and begin In 
column 1* The Input associated with some of the 
commands Is prompted for In conversational processing 
and In batch processing Is to be on cards Immediately 
following the command card. All commands except HALT 
may be Issued more 'than once. This Is especially 
useful In the conversational mode for correcting data 
entered In error. Many commands have opposite commands 
(I.e. STORE and NSTORE); the one given most recently 
Is the one that governs. The commands are given below 
In alphabetical order: 


COJltLNLIE 

EFFECT. ... Causes processing to continue In the batch 
mode Vi/hen an erroneous command Is 
encountered. Otherwise execution will 

terminate. In the conversational mode an 
Invalid command message will be Issued and 
processing will continue. 


PCM — £ielete j^hordwlse Modes 
Input NMAX2 

Effect. ,W1 1 1 limit the number of chordwise modes to 
NMAX2. NMAX2 must be 1 any previously 
entered NMAX2 and must be ^ the original 
number of chordwise modes/ NMAX, 
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DCR ^^Qel ete Xlhordwlse control £olnts 

I npu t • * • * • ( POEL) ““an Integer array of not more than 17 
points to be deleted. The array must be In 
ascending order. The indices referred to are 
always the original Indices (I.e, If there 
were 5 chordwise control points and the user 
deleted the 2nd and 4th^ then the Indices of 
the remaining points would be 1/ 3, 5 and not 
1/ 2, 3), 

Effect .... The chordwise control points specified are 
deleted. 


DPP — Helete a particular Point 

I npu t . . . . . P,M (P and M are I ntegers )“~These are 
always the original Indices. P denotes the 
chordwise control point number (P = 1 denotes 
the point closest to the leading edge). M 
denotes the spanwise control point number (M 
® 1 denotes the spanwise control point 

closest to the right side edge). 

Effect. ... The control point with the Indices P^M will 
be deleted. 


D-SM — Q,e\ete Spanwise Modes 

I nput. « • . • KK2; The new maximum number of spanwise 
modes. The number Includes all modes, not 
just the ones being used. For example If 
KK2 were 11, then this would Include 6 
symmetric and 5 antisymmetric modes, KK2 
must be to any previously entered KK2 
and/or ± KK, 

Ef feet , , . , L Iml ts the new maximum number of spanwise 
modes to KK2. 
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DSP -"De1 ete ^panwlse Control £olnts 

I npu t , . . . , (MDEL) --an Integer array of not more than 17 
points to be deleted* These always refer to 
the original control points and must be In 
ascending order. 

Ef feet .... The spanwlse control points with Indices 
In array (MDEL) will be deleted. 


HALT 

Ef feet .... Terminates execution. 


Effect. ... Begi ns a new wing leaving the following 
commands in effect: STORE or NSTORE/ SUM or 

NSUM, PCP or NPCP, RHS or NRHS, CONTINUE, 
PIND or NPIND. 


NEW 

Ef feet .... Begi ns a new wing resetting all commands and 
variables to their default settings. 


NPC 

Ef feet .... Reverses effect of PCP command; the control 
points will not be printed. (Sets logical 
variable POINTS to .FALSE,). This command 
need not be given unless the PCP command was 
previously given since this command reflects 
the default status. 


mm 

Ef feet .... Reverses effect of PINO command; Indices of 
rows and columns retained will not be 
printed. Sets logical variable PIND to 
.TRUE, This command need not be given unless 
the PIND command was prev i ous 1 y gl ven since 
this command reflects the default status. 
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iiPS 

Ef feet. Reverses the effect of PSOL command; the 
solutions will not be printed. This command 
need not be given unless the PSOL command was 
previously given since this command reflects 
the default status. 


MHi 

Ef feet Reverses the effect of RHS command; the right 
hand sides will not be printed. This command 
need not be given unless the RHS command was 
previously given since this command reflects 
the default status. 


K$TQ .RK 

Ef feet Reverses effect of STORE command““Sol u t 1 ons 
will not be stored on solution file. Since 
STORE Is a default this command must be 
issued to prevent solutions being stored. 


mm 

Ef feet Reverses effect of SUMMARY command--a summary 
of commands will not be printed. This 
command need not be given unless the SUM 
command was previously given since this 
command reflects the default status. 


£Q£»«£rint Control £oints 

Ef feet Control points will be printed. 
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PI ND --Prl nt INDlces 

Effect. .. .The Indices of rows and columns retained will 
be printed. 

PSOL — £rlnt SOL utl ons 

Effect. .. .Solutions will be printed. 


RCP — Restore Control joints 

Effect. .. .Restores the control points to that set 
originally on the Influence matrix and 
boundary condition files. This command can 
be used If the conversational user accidently 
specifies the wrong control points with the 
DCP or DSP commands. 


RHS — £lght Hand Jildes 

Effect. .. .Right hand sides will be printed. 


pPM — Restore pressure Modes 

Effect. .Restores the pressure modes. This command 
can be used If the conversational user 
accidently gives a- smaller value for NMAX2 or 
KK2 (DCM and DSM commands) than Intended. 


STAfll 

Ef feet ,,,. Causes processing to begin. This Is the last 
command that can be given for a particular 
solution. After the solution has been 
obtained and (If NSTORE has not been given) 
stored/ the user may delete more rows and/or 
columns and obtain more solutions. 
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mEL 

Effect. . 

lUMMEI 
Effect. . 

ISl 

Input... 
Effect. . 


.Solutions will be stored on solution file. 
This command need not be given unless the 
NSTORE command was previously given since 
this command reflects the default status. 


.A summary of commands will be printed. 


•A TSS command of 80 characters or less. 

•The command Is passed to the AMES* TSS 
operating system. After the system processes 
the command/ control returns to the program. 
This command Is a special one for the AMES* 
TSS version of the program. 
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7 SAMPLE CASE 

A sample conversational terminal session on the Ames* 
360/67 TSS computer system Is reproduced In this 
section with additional comments added In parenthesis. 
During this session 5 sets of solutions were 
determined. Each set consisted of 2 symmetric cases 
and 1 antisymmetric case. Each of these solution. sets 
was derived from the same Influence matrix to 
Illustrate how one can study the convergence of the 
kerhel function method In a very efficient manner (l.e. 
Influence matrices which are submatrices of larger 
Influence matrices are not calculated). 

For the first solution set the entire Influence matrix 
was used CI,e./ no control points or pressure modes 
were deleted). The chordwise control points used and 
the numbering of these control points are shown In fig, 
1 while the spanwise control points and numbering are 
shown In fig. 2. For this case the total number of 
control points equalled the total number of pressure 
modes. (30 for the symmetric cases and 25 for the 
antisymmetric case) 

For the second solution set the number of spanwise 
pressure modes was decreased from 11 to 5 (3 symmetric 
and 2 antisymmetric modes) while no control points were 
deleted. The solution thus obtained was a solution In 
the least squares sense (30 or 25 equations In either 
15 or 10 unknowns). 

For the third solution set the odd numbered spanwise 
control points were deleted leaving the spanwise 
control points shown In fig* 3, Since the effect of 
row and column deletion are accumulative/ this left as 
many control points as unknowns (15 for the symmetric 
cases and 10 for the antisymmetric cases). 

For the fourth solution set the second and fourth 
chordwise control points were deleted leaving chordwise 
control points on the leading edge/ mid-chord/ and 
trailing edge. Also the number of chordwise pressure 
modes was reduced from 5 to 3 leaving as many equations 
as unknowns (9 for the symmetric case and 6 for the 
unsymmetric case). 

For the fifth case It was desired to consider 3 
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chordwise control points and chordwise pressure modes 
and 11 spanwlse control points and spanwise pressure 
modes* For this It was necessary to restore all the 
control points and pressure modes using the RCP and RPM 
commands* Then the second and fourth chordwise control 
points were deleted and the number of chordwise 
pressure modes was reduced to 3/ yielding the desired 
matrix. 

The, output from the following session is given In 
Appendl x M ; 

LOGON user I d^ password/ termi nal Id 

AMES USYSLIB 

JOBLIBS SYSULIB 

JBLB MEDAN 

CALL S0LVR$ 

ENTER BATCH 
(carriage return) 

ENTER ODISK 
1 

OUTPUT IS ON .**0UTPUT,S0LVR.N1... 

CANCELLED! DDNAME FT04F001 UNKNOWN 
(Messages such as the above occur because of the 
automatic file defining feature and do not Indicate any 
error , ) 

ENTER IDDl/ I DD2/ IDD3 
4 2 9 

CANCELLED! DDNAME FTllFOOl UNKNOWN 
CANCELLED! DDNAME FT08F001 UNKNOWN 
ENTER VALUE FOR NTYPE/ 0 FOR SYMMETRIC, -1 FOR 
ANTI “SYMMETRIC, 1 FOR BOTH 
1 
• 5 * 

PCP 

PSOL 

RHS 

+ 

SUMMARY 

+ 

TSS 

CPUTIME? 

4.350 SECONDS 

START 
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ID4= 15 

CANCELLED: DDNAME FT12F001 UNKNOWN 

♦ 

TSS 

CPUTIME? 

9.208 SECONDS 


DSM 

ENTER NEW MAXIMUM NUMBER OF SPANWISE MODES 
IS 11 

ORIGINAL VALUE WAS 11 
5 


PIND 

+ 

START 
|D4» 16 

+ 


CPUTIME? 

11.996 SECONDS 

+ 

DSP 

ENTER SPANWISE CONTROL POINTS TO BE DELETED 
ASCENDING ORDER 
13 5 

+ 

START 
ID4» 17 
+ 

CPUTIME? 

14.072 SECONDS 

+ 

DCP 

ENTER CHORDWISE POINTS TO BE DELETED 17 
ASCENDING ORDER 

2 4 

+ 

DCM 

ENTER NEW NUMBER OF CHORDWISE MODES OLD 
ORIGINAL VALUE VMS 5 

3 

+ 

START 
ID4» 18 
♦ 


OLD MAXIMUM 


17 MAXIMUM 


MAXIMUM IN 


VALUE IS 5 
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TSS 

CPUTIME? 

15,956 SECONDS 

+ 

RCP 

+ 

RPM 

+ 

DCP 

ENTER CHORDWISE POINTS TO BE DELETED 17 MAXIMUM IN 
ASCENDING ORDER 
2 4 


DCM 

ENTER NEW NUMBER OF CHORDWISE MODES OLD VALUE IS 5 
ORIGINAL VALUE WAS 5 
3 

+ 

START 
ID4« 19 

TSS 

CPUTIME? 

18.529 SECONDS 

+ 

HALT 

PROGRAM TERMINATED BY HALT COMMAND 
TERMINATED: STOP 

PRINT OUTPUT. SOLVR.Nl, PRTSP»ED I T^ STATI ON-RMT05 
PRINT BSN-2232/ 300 LINES 

LOGOFF 
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8 SUBROUTINES 


8.1 STOP2 This subroutine prints an error message and 
terminates execution. It assumes an 80 
character message with words of 4 characters 
each. If the message from the calling 
program Is not this long/ then some garbage 
will appear at the end of the message. This 
program will require revision In the FORMAT 
statement If the program Is used on a system 
with word length differing from 4 characters. 


8.2 SOLVIT Equation solving subroutine which uses 
Householder's method described In section 3. 
Also see the program listing for more 
documentation. 


8.3 KRUNCH 


This subroutine Is 
arrays ( I NDCOL) and 
Indices of rows and 
See the program 
expl anat I on. 


used In determining the 
(INDROW) which are the 
columns to be retained. 
1 Isting for further 


8.4 AIMFIL 


This subroutine Is used In the AMES' TSS 
version only. Its purpose Is to Issue a 
RELEASE command and OOEF command for the 
Influence matrix file. 


8.5 BCFIL This subroutine Is used In the AMES' TSS 
version only. Its purpose Is to Issue a 
RELEASE command and a DDEF command for the 
boundary condition file. 


8.6 BNKFIL This subroutine Is used In the AMES' TSS 
version only. Its purpose Is to Issue a 
RELEASE command and a DDEF command for the 
solution file. 
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8.7 OBEY 


8.8 CVRT 


This Is an Ames* library subroutine used for 
passing commands from the program to the TSS 
operating system. 


This Is an Ames* library program used for 
preparing RELEASE and DDEF commands for use 
In OBEY. 
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9 LNPJIT ElLEl 


The following disk files are read by the program. 
For a complete description of the Information on 
the files see refs. 2-3. The AMES* TSS version of 
the program issues Its own DDEF commands for the 
flles/ so none need be given. For other systems 
appropriate control cards will have to be supplied 
for units 8/ 9^ 11/ and 12. 


9,1 Aerodynamic Influence Matrix File (AIM file) 


This file Is a variable record length file and is 
read from unit 11. 


The first record contains Identification and title 
Information plus Information about the size of the 
matrix and location of spanwlse and chordv;lse 
control points. 


The second and subsequent records contain the 
Influence matrix Itself. This file Is generated 
by the Influence matrix program (ref. 2), 

On the Ames* TSS system this file has the name 
AIM.XI.XJ where I Is the numerical value of IDl 
and J Is the numerical value of ID2. 


9.2 Boundary Condition File (BC file) 


This file Is a variable record length file read 
from unit 8. 


The first record contains Identification and title 
Information plus Information Identifying the type 
and number of symmetric and antisymmetric cases. 


The next NSYM records are right-hand sides for 
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symmetric cases where NSYM equals the number of 
symmetric cases. The next NASYM records are the 
right-hand sides for antisymmetric cases where 
NASYM equals the number of antisymmetric cases. 
In the case of an unsymmetric wing there will be 
NSYM + NASYM right-hand sides. 

This file Is generated either by the geometry 
program (ref. 2) or the boundary condition program 
(ref. 4). 

On the Ames* TSS system this file has the name 
BC.XI.XK where I Is the numerical value of IDl and 
K Is the numerical value of ID3/ which Is 
determined from IDFILE at the time the program Is 
run and Is found In the program output. 


9.3 Identification Number File 


This file Is read from unit 9 and rewritten on 
unit 9. 


This file contains Identification mumbers In 
binary form. The fourth number of the file Is the 
one used to Identify the solution. Whenever 
solutions are to be stored the file Is read and 
updated. This new Identification number Is then 
written on the solution file. 

On the Ames* TSS system this file has the name 
IDFILE. 
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MATRIX MANIPULATION AND EQUATION SOLVING PROGRAM 


10 ^UTLON £LL£ 


This Is a binary file with variable record lengths and 
Is written on unit 12. The AMES* TSS version Issues 
Its own DDEF commands (control cards) for the file so 
none need be given by the user. On the Ames’ TSS 
system this file has the name BNK.XI .XJ.XK.XL where J 
Is the numerical value of IDl/ J Is the numerical value 
of ID2^ K Is the numerical value of 103/ and L Is the 
numerical value of JD4. On another system appropriate 
control cards will have to be supplied for unit 12. 


FIRST Record 

IDl/ ID2/ ID3/ 104/ TITLE/ SYM/ NSYM/ (BCS( I )/ I -1/ 10 ) / 
(BCAS(I)/ 1=1/10)/ NMAX2/ KK2/ JJ/ PPNEW/ MMNEW/ NROWSA/ 
DELTAO/ EPS/ MACH 


IDl 

Identification number of the geometry 
file from which AIM file and BC files 
were derived. 

ID2 

Identification number from AIM file. 

ID3 

Identification number from BC file. 

ID4 

Identification number assigned to this 
solution. 

(TITLE(26)) 

Title Information. 

SYM 

Logical variable whose value Is .TRUE, 
for a symmetric wing and .FALSE, for 
an unsymmetrlc wing. 

NSYM 

Number of symmetric cases. 

NASYM 

Number of antisymmetric cases. 

(BCS) 

Ten element logical array Identifying 

the types of the symmetric cases. 

(BCAS) 

Ten element logical array Identifying 
the types of the antisymmetric cases. 
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MATRIX MANIPULATION AND EQUATION SOLVING PROGRAM 


NMAX2 

KK2 

JJ 


PPNEW 
MM NEW 
NROWSA 


DELTAO 

EPS 

MACH 


New number of chordwlse modes. 

New number of spanwise modes. 

The number of spanwise Integration 
points used by the Influence matrix 
program. This number comes from the AIM 
file. 

New number of chordwlse points. 

New number of spanwise points. 

The number of rows after deletion. 
This number Is equal to the number of 
control points actually used. If both 

symmetric and antisymmetric solutions 

have been obtained/ then NROWSA Is equal 
to the number of control points used for 
the symmetric cases. 

Parameter In the calculation of the 

1 nf 1 uence matr 1 x. 

Parameter In the calculation of the 

I nfl uence matr I x. 

Mach number (floating point). 


ji£ii Mm Rmmi 

(X(J)/J»1/ NROWSA) The nonzero coefficients/ Bnk/ of 

the expansion for the pressure for 
the symmetric cases. These are 
stored so as to be compatible with 
the following READ statement: 

READ(12)((BNK(N/K)/N»1/NMAX2)/K*1/ 

KK2/KJUMP) 

where KJUMP»2 If the wing Is 
symmetric or XJUMP^l If the wing Is 
unsymmetr 1 c. 
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MATRIX MANIPULATION AND EQUATION SOLVING PROGRAM 


mi mm 

(X( J)/ J«l/ NROWSA) The nonzero coefficients/ Bnk/ of 

the expansion for the pressure for 
the antisymmetric cases. These are 
stored so as to be compatible with 
the following READ statement: 

READ(12)((BNK(N/K)/N=1/NMAX2)/K- 

KJUMP/KK2/KJUMP). 
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MATRIX MANIPULATION AND EQUATION SOLVING PROGRAM 


11 ££££ R£ N, PES 


1 . 


Medan^ R. T. : 
Method of 
Theory. 
1974. 


Improvements to the Kernel Function 
Steady^ Subsonic Lifting Surface 
NASA Rept. No. TMX-62,327, Mar. 


2 . 


MedaO/ R. T.: Geometry Program for 

Lifting Surface Theory. NASA 
TMX-62y309, Sept. 1973. 


Aerodynamic 
Rept. No. 


3. Medap/ R. T./ 
Program 
Theory. 
1973. 


and Ray/ K, S. ; 
for Aerodynmalc 
NASA Rept. No. 


Influence Matrix 
Lifting Surface 
TMX 62/324/ Dec. 


4. Medan/ R. T./ and Ray, K. S. : Boundary Condition 

Program for Aerodynamic Lifting Surface 
Theory/ NASA Rept. No. TMX-62/323/ Dec. 1973. 


5. Medan/ R, T,/ and Ray/ K. S,; Normal Loads Program 
for Aerodynamic Lifting Surface Theory/ NASA 
Rept. No. TMX-62/326/ Feb. 1974. 


6. Medan/ R. T./ and Ray, K. S. : Plotting Program 

for Aerodynamic Lifting Surface Theory/ NASA 
Rept. No, TMX-62/321/ Nov. 1973. 


Tullnlus/J.; Clever/ W.; Niemann/ A.; Dunn/ K.; 
and Gaither/ B.: Theoretical Prediction of 

Airplane Stability Derivatives at Subcritical 
Speeds. Rept. no* NA-72-803/ North American 
Rockwell Corp./ 1972. 
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FIGURES 







Figure 2 - The Spanwlse Control Point Distribution and 
Numbering Scheme for the Sample Case. (Note: The 
wing In the figure Is not the wing used In the sample 
case. ) 
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OUTPUT FROM SAMPLE CASE 
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DETER.HINAT ION OF BNK COEFFICIENTS 
PFCTANGULAR WING AR - ? 11-13-73 


101 

= 

4 

102 

s 

2 

103 

s 

9 

104 

3 

15 

JJ 

ac 

191 

UNSYM 

3 

0 

NODES 

3 

1 

ntype 

SS 

1 

PPNEW 

3 

5 

^MNFW 

3 

11 

NMAX2 

. = 

5 

KK2 

s 

11 

OELTAO 

3 

4.0000 

MACH 

X 

0.0000 

PCS 

3 

TTFFFFFFFl 

BCAS 

3 

FTFFFFFFP 

SUMMARY 

OF 

COMMANDS 


PCP 

PSOL 

PHS 

SUMMARY 

TSS 

START 

CHOROWISE CONTROL POINTS USED 

€♦00000 
O, 14645 
0.50000 
0.85355 
1.00000 


SPAMWISE CONTROL POINTS USED 


0.96593 

0.06603 

0.70711 

0.50000 

0.25082 

- 0.00000 

BOUNDARY CONDITIONS AT ORIGINAL CONTROL 
(L.E. TO T*e.. RIGHT TIP TO LEFT TIP) 


CASE I 


1.00000 

1.00000 

1. 000 00 

1.00000 

1.00000 

1.00000 

1. DOOOO 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

l.OOOOC 

1.00000 

1.00000 


POINTS 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 


1. 00000 


1.00000 



CASE 7 


0.00000 

0.14645 

0,50000 

■ 0.85355 

1.00000 

0 . coooo 

0.14645 

0.50000 

0.85355 

1.00000 

0.00000 

0.14645 

0.50000 

0.95355 

1.00000 

o.ooooc 

0.14645 

0.50000 

0.85355 

1.00000 

0.00000 

0.14645 

0.50000 

0.85355 

1. COOOO 

0.00000 

0.14645 

0.50000 

0.85355 

1.00000 

NUMBER Of 

control POINTS * 

30 



NUMBER OF 

PRESSURE MODES = 

30 




THE DETERMINANT OF ( A-TRANSPOS E )*A * 0.11T73389E 22 


SOLUTION FOR SYMMETRIC CASES 


CASE 1 

1 0.9592875 

2 -0.3A34652 

3 -0*0878328 

4 -0.0263244 

5 -0.0096903 


CASE 2 

0.29 37329 
0.6579308 
-0.0227847 
0.0056143 
-0.0019472 



6 0.1084463 0.0032115 

7 -0.1937814 0.0297757 

8 -0.1148094 -0.0201028 

9 -0.0630316 0.0021874 

10 -0.0266566 -0.0050513 

11 0.0439750 0.0042856 

12 -0.0883476 -0.0059570 

13 -0.0770952 -0.0120187 

14 -0.0612980 -0.0035019 

15-0.0312412 -0.0055371 

16 0.0233391 0.0031349 

17 -0.0468022 -0.0057873 

18 -0.0452612 -0.0070768 

19-0.0428134 -0.C045091 

20 -0.0246244 -0.0042022 

21 0.0119725 0.0017669 

22 -0.0239923 -0.0034023 

23 -0.0237498 -0.0037836 

24-0.0236830 -0.0030724 

25 -0.0144601 -0.0024254 

26 0.0036820 0.0005621 

27 -0.0073745 -0.0010958 

28 -0.0073520 -0.0011825 

29 -0.0074010 -0.0010369 

30 -0.0046311 -0.0007719 

BOUNDARY CONPITTHN^ AT ORIGINAL CONTROL POINTS 
(L.E. TO T.E.. RIGHT TIP TO LEFT TIP) 


0.96593 

0.96593 

0.96593 

. 0.96593 

0.96593 

0.86603 

0.86603 

0.86603 

0.86603 

0.86603 

0.70711 

0.70711 

0.70711 

0.70711 

0.70711 

0.50000 

0.50000 

0.50000 

0.50000 

0.50000 

0.2588? 

0.25882 

0.25882 

0.25882 

0.25882 

NUMBER OF 

CONTROL POINTS = 

25 



NUMBER OF 

PRESSURE MODES = 

25 




THE OETERHfNANT OF ( A-TRANSPOSE A = 0.384417?3E 18 


SOLUTION FOP ANTI-SYMMETRIC CASES 


CASE 1 

1 0.3958588 

2 -0.3087006 

3 -0.1107676 

4 -0.0433692 

5 -0.0158296 

6 0.C606708 

7 -0.1169108 

8 -0.0854544 

9 -0.0566794 
10 -0.0255904 



11 0,027072l 

12 -0.0544521 

13 -0.0503660 

14 -0.3443864 

15 -0.0238526 

16 0.0141472 

17 -0.0283824 

18 -0.0277684 

19 -0.0272557 

20 -0.0161263 

21 0.0062318 

22 -0.0124904 

23 -0.0124014 

24 -0.0124801 

25 -0.0077C80 

SOLUTION FOUND A NO STORED 


FINISHED 


Df:Ti:RMlNATIOM -JP BNK COEPFICIFNTS 
RECTANGULAR WING AR * 2 11-13-73 


U3 


ini 

102 

ID3 

104 

JJ 

UNSYM 
MODES 
NTYPE 
PPNEW 
MMNEW 
MM AX 2 
KK2 

OEL^AO 

MACH 

BCS 

BCAS 


4 
2 
9 

16 

191 

0 

1 

1 

5 
II 

5 

5 

4.0000 

0.0000 

TTFETFFFFF 

FTPFFFFFFF 


SUMMARY OF COMMANDS 

TSS 

PSM 

PINO 

start 

CHORDWISE CONTROL POINTS USED 

0.00000 
0.14645 
0.50000 
C .85355 
1. 00000 



SPANWISE CONTROL POINTS USED 

0.96593 

0.06603 

0.70711 

0.50000 

0.25882 

- 0.00000 

INDICES OF COLUMNS TO BE RETAINED 


1 

2 

3 

4 

5 

13 

14 

15 

21 

22 

25 





INDICES 

OF 

ROWS TO 

BE RETAINED 

1 

2 

3 

4 

5 

8 

9 

10 

11 

12 

15 

16 

17 

18 

19 

22 

23 

24 

25 

26 

29 

30 




NUMBER 

OF 

CONTROL 

POINTS = 

30 

NUMBER 

OF 

PRESSURE 

MOOES = 

15 


THE DETERMINANT OF < A-TRANSPOSE 1*A 


SOLUTION FOR SYMMETRIC CASES 


ll 12 

23 24 


6 7 

13 14 

20 21 

27 20 


0.92918963E 09 


CASE 1 

CASE 2 

1 C. 0563015 

0.2933132 

2 -0.3391642 

0.6563294 

3 -0.085B390 

-0.0223944 

4 -0.0247415 

0.0057314 

5 -0.0089764 

-0.0018154 

6 0.1117279 

0.0036379 

7 -0.2003296 

0.0290937 

e -0.1217275 

-0.0213837 

9 -0.0688910 

0.0018656 

10 -0.0288136 

-0.0055897 

11 0.0466561 

0.0046587 

12 -0.0936345 

-0.0063517 

13 -0.0796093 

-0.0128674 

14 -0.0597491 

-0.0029139 

15 -0.0285094 

-0.0053198 


INDICES OF COLUMNS TO BE PETAINFO 


6 

T' 

8 

9 

10 

16 

17 

18 

19 

20 





INDICES 

OF 

ROWS TO 

8E RETAINED 



1 

2 

3 

4 

5 

6 

7 

8 

Q 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 




NUMBER 

OF 

CONTROL 

POINTS = 

25 



number 

OF 

PRESSURE 

MOOES = 

10 




THE OETERMINAMT OF < A-TR ANSPOS E ) ♦ A = 


0.74841275E 06 



SOLUTION FOR ANT I -SYMMETRI C CASES 


CASE 1 

1 0,3965770 

2 -0.3107769 

3 -0.1147077 

4 -0.0454141 

5 -0.0156609 

6 0.066P647 

7 -0.1312739 

8 -0.0916794 

9 -0.0548593 
10 -0.0220116 

SOLUTION FOUND AND STORED 


FINISHED 



determination of rnk coefficients 

RFCTANCULAR WING AP * 2 11-13-73 


lOl 

= 

4 

102 

5 

2 

103 

S 

9 

ID4 

= 

17 

JJ 

S 

191 

UNSYM 

= 

0 

MOOES 

= 

1 

NTYPE 

s 

1 

ppnfw 

= 

5 

MMNEW 

s 

5 

NMAX2 

s 

5 

KK2 

= 

5 

DELT AO 

= 

4.0000 

MACH 

3B 

0.0000 

PCS 

* 

ttffffffff 

BCA5 

s 

PTFFFFFFFF 

SUMMARY 

OF 

COMMANDS 


TSS 

DSP 

START 

CHQROWISE CONTROL POINTS USED 

0,00000 

0,14645 

0.50000 

0.05355 

l.OOOCO 



SPANWISE CONTROL POINTS USED 


0.86603 

0.50000 

- 0.00000 



INDICES 

OF 

COLUMNS 

TO 

BE 

RETAINED 




1 

2 

3 


A 

5 

ll 

12 


13 

lA 

15 


21 

22 

23 

2A 


25 









INDICES 

OF 

ROWS TO 

BE 

RET 

AINEO 



-p- 

6 

7 

8 


o 

10 

16 

17 


18 

19 

20 


26 

27 

28 

29 


30 









NUMBER 

OF 

CONTROL 

POINTS 

* 15 




NUMBER 

OF 

PRESSURE 

MOOES 

* 15 




THE DETERMINANT OF ( A-TRANSPO$ E A = 0.1782394CE 07 


SOLUTION FOR S 

YMMETRIC CASES 

CASE 1 

CASE 2 

1 0.95A962A 

0.2Q30A9 8 

2 -0.33A5991 

0.6591860 

3 -0.0790099 

-0.0212807 

A -0.0175583 

0.0067213 

5 -O.OOA7A1P 

-0.0010586 



6 0.0948424 0.0012067 

7 -0.1665186 0.C3 35271 

8 -0.0860C90 -O.C 156435 

9 -0.0367432 0.0053368 

10 -0.0117176 -0.0023978 

11 0.0195498 0.0010045 

12 -0.0393625 0.0000173 

13 -0.0300078 -0.0045080 

14 -0.0172923 0.0009768 

15 -0.0065783 -0.0012288 

INDICES OF COLUMNS TO BE PETAINED 

67 8 9 10 

18 19 20 

■p* 

vn 

INDICES OF ROWS TO BE RETAINED 

6 7 8 9 10 

18 19 20 

. NUMBER OF CONTROL POINTS * 10 

NUMBER OF PRESSURE MODES = 10 


THE DETERMINANT OF ( A-TRANSPOSE )*A 


SOLUTION FOP ant I-SYMMETP.IC CASES 


CASE 1 

1 0.3885300 

2 -0.2940048 

3 -0.0962846 

4 -0.0290404 

5 -0.0077700 


16 17 


16 17 


0.86715039E 04 



6 0.0452002 

7 - 0.0878599 

8 - 0.0553636 

9 - 0.0276744 
10 - 0.0093136 

SOLUTION FOUND AND STORED 


FINISHED 



PETFRMINATION GF BHK COEFFICIENTS 


RfcCTANGULAR WING AR 

lOl 

s 

4 

102 

= 

2 

103 

s 

9 

104 

s 

IR 

JJ 

= 

191 

UNSYM 

s; 

0 

MUDES 

s 

1 

NT YOF 

= 

1 

PPNEW 

= 

3 

MMNFW 

s 

5 

NMAX2 

s 

3 

KK2 

« 

5 

OELTAO 

s 

4.0000 

MACH 

= 

0.0000 

BCS 


TTFFFFFFFF 

BCAS 


FT^FFFFFFF 


SUMMARY Of COMMANDS 


TSS 

OCR 

DCM 

START 

CHOF.DWISE CONTROL POINTS USED 

o.ooooc 

0*50003 
I • 00030 



SPANWlSe CONTFGL POINTS USED 

0.86603 

0.50000 

- 0.00000 

INDICES OF COLUMNS TO BE PETAINEO 

1 2 3 11 12 13 21 

22 23 

INDICES OF FOWS TO BE RETAINED 

18 20 26 

9 
9 


THE OETEPMINANT OF ( A-TRANSPOSE )»A * 0.10721305E 05 


6 0 10 16 

28 30 

^ NUMBER OF CONTROL POINTS = 

«> NUMBER OF PRESSURE MOOES = 


SOLUTION FOR SYMMETRIC CASES 


CASE I CASE 2 


1 0.9595866 

2 -0.3473339 

3 -0.0799000 


0.2921830 

0.6642576 

-0.0210B88 


4 0.1041752 

5 -O.10B5B46 

6 -0.0907550 


0.0008584 

0.0355538 

-0.0151825 


7 0.0233645 

8 -0.04^?085 

9 -0.0314681 


0.0011664 

-0.00017^7 

-0.0043278 



INDICES IF COLUMNS TO BE RETAINED 

6 7 8 16 IT 18 

INDICES OF ROWS TO BE RETAINED 

6 a 10 16 la 20 

NUMBER OF CONTROL POINTS =6 
NUMBER OF PRESSURE MOOES = 6 


THE determinant OF I A-TPANSPOSE )^A = 0.41941284E 03 


SOLUTION for ANTl-SYMMetRlC CASES 


CASE I 

1 0.3961701 

2 -0.3133379 

3 -0.098U25 

U 0.0517030 

5 -0.1026420 

6 -0.0576963 

SOLUTION FOUND AND STORED 


FINISHFD 



DE'^'ERM INAT ION OF BNK COEFFICIENTS 
RECTANGULAR WING AR = 2 11-13-73 

IDl » ^ 

102 * 2 

103 = ? 

10 ^ * 19 

JJ = 191 

UNSYM * 0 

MOOtS = 1 

NTYPE = 1 

PPNFW = 3 

MMNFW « 11 

N'MAX? = 3 

KK2 = 11 

DELTAO = 4.0000 

MACH = 0.0000 

RCS * TTFFFFFFFF 

BCAS = FTFFFFFFFF 


SUMMARY OF COMMANDS 

TSS 

PCP 

RPM 

DCP 

DCM 

START 

CHOROWISE CONTROL POINTS USED 

n. 00000 

0.50000 

1.00000 



SPANWISE CONTPGL PHINTS USED 

0,96593 

0,86603 

0,70711 

0,50000 

0.25882 

- 0.00000 

INDICES OE CHLUMNS TO BE RETAINED 

1 2 3 11 12 

22 23 31 32 33 

43 51 52 53 

INDICES OF ROWS TO RE RETAINED 

1 3 5 6 0 

13 15 16 18 20 

25 26 28 30 

NUMBER OF CONTROL POINTS * IP 
NUMBER OF PRESSURE MOOES = 18 


THE DETERMINANT GF ( A-TRANSPOSE )*A 


SOLUTION FOR SYMMETRIC CASES 


CASF 1 CASE 2 

1 0.9661451 0.2932169 

2 -0.3604457 0.6623058 

3 -0.C8B7532 -0.0225372 


13 21 

41 42 


10 11 

21 23 


0.62640890E 14 



4 

0.1228009 

0.0036634 

5 

-0.2260956 

0.D3G1290 

6 

-0.1176664 

-0.0194774 

7 

0.054r528 

0.0054766 

8 

-0.1107294 

-0.0082956 

9 

-0.0782859 

-0.0114408 

10 

0.0287693 

0.0039528 

11 

-0.0575409 

-0.0074988 

12 

-0.0435586 

-0.0C65170 

13 

0.0140763 

0.0021257 

14 

-0.0280615 

-0.0041540 

15 

-0.0213801 

-0.0033279 

16 

0.0041469 

0.0006469 

17 

-0.0082601 

-0.0012729 

18 

-0.0063017 

-0.0010010 


INDICES 

OF 

CtiLUMNS 

TO 

BE 

RETAINED 



6 

7 

8 


16 


17 

18 

26 

27 

29 

36 


37 


38 

46 

47 

48 









INDICES 

OF 

ROWS TO 

BE 

RETAINED 



1 

3 

5 


6 


B 

10 

11 

13 

15 

16 


18 


20 

21 

23 

25 

NUMBER 

OP 

CONTROL 

PC IN 

TS 

S 

15 



NUMBER 

OF 

PRESSURE 

MODES 

7 

15 




THE D6T6RMINAMT OF ( A-TRAN5P0SB )*A 


0.431 1908AF 12 



SOLUTION FOR ANT I-SYMMFTRIC CASES 


CASE 1 

1 0.4065139 

2 - 0.3341637 

3 - 0.1126727 

4 0.0720457 

5 - 0.1435348 

6 - 0.0877573 

7 0.0336897 

8 - 0.0677919 

9 - 0.0501169 

10 0.0170100 

11 - 0.0339469 

12 - 0.0258488 

13 0.0071492 

14 - 0.0142447 

15 - 0.0108606 

SOLUTION FOUND AND STORED 


FINISHED 



APPENDIX II 


COMPUTER PROGRAM LISTING 



iNtEOE« UCl,U(»,UcO,U6,U20#Un ^UU^UNSYMjPPiPiC^^TYPg^SWTVPE, 
ipd€i,hold»pp^ew,odisk 

IOGICaL BCS.BCASf INVEPS^RtAOA, B atch, store, RHS,C 0MSU«, 
lPOlNTS,cONTIN,cOHV,PSOU,PINO,SKlP,cHtCK,OPEN^,SYM, aSYM, 
2BCS2,BCASa 
PEAL Mach 

PEaL HaL$,NEW$,KNE$,NST$,NSU$,NPc$,NRH$,NPSS,NP1$ 

OIME^^SION PDEL(17) jHnEL(l7),TITLE(26)f INOPOW(IOO)# 

1 INOcOLf iOO)f 1 NOOEUHOO)#CaPOIMc«^ 100) ^ETACP(ar) I 
2NlNDEXC^T),CHlCP(20)»AlFA(«0,l0)#BCSfl0),BCAS(10># 

SXET ACP(P7),A( 1000), A«< 120)# 

^11(120)# At4(B0#P0)#XCHlCP(20),RM8RE8(20) #X(60# 10) # 

SAlNV(l# 100 )#Ll 5 T 3 f 2 # 25 )#Ll 8 Tl( 20 )#Ll 5 T 2 («r)fBCs 2 ( 10 )#dCAs 2 f 10 ) 


C*».*.T8ScOM is OHLY used TO STORE AWESt TS8 COMHaNDS, 

C 

OIMEHSION TSScOM(20) 

C 

equivalence (A1NV(1, 1)#XETaCPc 1))# (MOEL(l)#PDEL(l))f 
l(XCHICP(l)#Afl)) 


C 

C 

C 

C 

C 

C 

C 

C 

C 


♦•..♦NX SHOULD hE the first DIMENSION OF (X) aNO SHOULD pE 

large as the number of MOOES TOTaL, the SECOND OIH. 
,,,,,ENSiON Of (X) SmOUlD BE a8 LARgF aS ThE NUMBER OF CaSES 
•••••(BOTH SYMMETRIC AND ANTISYMMETRIC), 

NB SHOULD rE THE FIRST DIMENSION OF (ALFA) AND SHOULD pE 

....•AS LARg^ as ThE number of CONTROL POtNTS BEFORE ANY 
,,,,,aRE deleted, the second dimension of (ALFA) SHOULD BE AS 
,,...LAHGE as THE number OF CaSES BEING TREaTEO, 

DATA UCl/S/#UCO/6/, US/S/, U20/OR/,Ull/ll/#Ul2/l2 /, blank /1h 
ICASE/RHCaSE/, IFMT/0/# NX/80/, NB/SO/, MAIN V/1/, 

2lNVER5/,FALSE./#REAOA/.TRUE./iNcOMK/ioo/ 

DATA HALS / ShHAL / 
data news / 3HNEW / 

DATA FNE$ / IJHKNE / 

DATA OCMS / JhDCM / 
data DSM$ / 3HDSM / 
data DCPS / 5HDCP / 




/# 
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data 0SP$ / 3H0SP / 

data 0PP$ / 3H0PP / 

data cons / 3HCQN / 

Data stas / shsTa / 

OATA NSTS / 5MN8T / 

OaTa 8T0S / 5HST0 / 

OAtA SUNS / 3 HSUh / 

data NSUS / 3HNSU / 

OATA PCPS / 3HPCP / 

DaTa NPCS / 3HNPC / 

OATA PHSs / 3HRHS / 

data NRHS / 5HN«H / 

data PSOS / 3MPS0 / 

data NPSS / 3HNPS / 

data pins / 3HPXn / 

OATA NPls / 3HNP1 / 

data StTs / 3HSET / 

data TSSS / 3HTS8 / 

OaTa «PmS / 3HPPm / 

OATA HCPS / 3NWCP / 

C 

OPrNA a .PAtSI, 

COM 8 blank 
MRIT e(UfOrAOl) 

C 

C establish WHETHER THIS IS A BATCH OR CONVERSATIONAL JOB 

READrUCI#799) BATCH 

conv • , NOT, batch 

C 

j7 CONTINUE 
I NO « fj 

IND t 1 WHEN BOTH symmetric AND ANTISYMMETRIC SOLUTIONS 

are to be obtained 

a when only one type OF SOLUTION IS TO bE 

obtained, 

LISTI(I) » 0 
L1ST2(1) • 0 
ASYM a , false, 

IF(OPEN^») ENpiFILE U6 
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IFrOPENb) CALL OBEvnbi UH«EL6 aSE FTOaPOOJ ) 

U6 9 ucn 

IF(CONV) WRn£(UrO>8|5) 

«EAOrUCI#8lO)OOISK 
IE(001SK,NE,0) U6 a 
lFfU6,NE,UcO^ OPEMfea.TRUE, 

Iff . nOt.OPcN^) go to 18 

aMES» TSS system only, Ob6y SENDS COMMANDS TO TSS 
• •,.,CVRT IS A library program FOR IN-CORE RE AO , RE w 1 NO , NR 1 TE , 
OOlSK a MOD(MAK0n,OoiSK), tO) 

WRITE( UCO,RSn ODI8K 

Call obey 1 1<>, ishrelease fto4ifooi ) 

CALI CVRT(OolSK, I# 

I a«H(roOEF FTo«F001#»OUTPUT,SOLVR,Ni,nibiO f 

call 0BEY(32,X) 
rewind a 
l8 CONTINUE 

ENTER 10 numbers 
IF fCONv) wRlTtfUCOiBEO) 

READfUCNSlo^ lODl # 1002# ID03 
lF(IODl,Nt“,0} CHgCK a ,TRU^, 

CHECK IF False suppresses checking of aim and 
bc file identification numbers 

SKIP B , FALSE, 

SKIP when ,TRUE. SUPPRESSES READING OF BOUNDARY CONDITION 
FILE— IT IS SET TO TRUE WHEN A CASE IS COMPLETED AND 
MORE SOLUTIONS aRE TO BE OBTAINED BY DELETING ROWS AND 
COLUMNS 

OPEN AND READ AIM FILE 

ON standard computer remove FOLLOWING CALL AND OPEN 
statements and ATTACH STATEMENT NUMgER 2o TO REaD 
statement 
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,FOW aMES» TSS version only, AIMFIL issues DOEf TO OPEN 
.AIM file, 

CALL AIMFlLdOOl# 1002) 

RtAO(Un) I0ldO2,NSEQ.TlTLE,NT!TL#PP#rWTYPE,MM,MREF, 
lSi>.TYPe,NMAX,KK,HOOES,UNSYM,LMlN,LMAX,JJMAX, JJ.MaCH.EPS, 
2DELTA0dCHICP(ndpl#^P)dNlN0exM)f 
5(ETACP(M),M«1,MM) 

PP • number op CHOROwISE control points 

MM • number of SPANWISE CONTROL POINTS ON ENTIRE MING 

MOOES • Types of modes available 

0 B symmetric only 

less Than O s ANTfSYHMETRjC ONLY 
MORE THAN 0 » ALL MODES AVAILABLE 
UNSYM » SYMMETRY INDICATOR 

NOT 0 « UNSVMMETRIC WING 
f) « SYMMETRIC WlNG 

CMICP « LOCATIONS OF CHOROWISE CONTROL POINTS 
ETaCP • LOCATION OF 5PANWISE CONTROL POINTS 

ASYM m .true, FOR UNSYMMETRlC WING 

• .false, for SYMMETRIC WiNG 
SYm e .TRUE, FOR SYMMETRIC wjNr, 

• .false, for UNSYMMETRlC wlNG 

IF (UNSYM, NE ,0) ASYM » ,TRUE, 

SYM a .NOT, ASYM 

unless Check ■ .false. ChEck to be sure id on Files 
Match those rEaO in 

.THIS checking is not NECESSaRY FOR THE aMESi VERSION. rUT IS 
.DONE anyway, 

IF ( .NOT. CHECK) gO TO 
IF(Ini.EQ.IDDl) GO TO 

IF (BATCH) Call STOP;>(UCO. f iDDi DOES NOT MaTCH aIM FlLEw 
IFLOaT(IOD) 

pause » lODl DOES NOT MATCH AIM FILE VaLUEi 
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55 lFfI02.tO»ID02) GO TO 45 
IF (BATCH)CAUI 

IST0P2CUC0»^ 100? REiO IN DOES NOT MATCH AIM VALUE W 
2FL0ATC10?)) 

PAUSE • lOD? DOES NOT MATCH AIM FILE VaLUE ' 

C 

C OPEN boundary condition FILE 

45 CONTINUE 

.....BCFIL ISSUES OOEF CONTROL CaRD FOR AMES* T8S VERSION, 
Call -BCFILflOtf IDD5) 

RE AO (Us) 10UlD3rTlTLE,UNSYM,N8YM#NASYM, (BCS(I)# 

Half 10), (BCAS(I),l«l#lO)#PP#CWTYPE,MM,MMP,NHPA,MReF, 
25HTYPE, (XCHlCP(I),IalfPP)f (NlNDEX(I)rI»UMMP), 
5(XETACP(I),Ial,MMP^ 

N$YM 9 NUMBER OF SYMMETRIC CASES 

nasyh • number of anti-symmetric CaSES 

MmP c number Of SPANNSIE CONTROL POINTS USED fOR 

THE symmetric Cases 

MMPA » number of SPaNWISE CONTROL POINTS USED FOR 
THE anti-svmmetric Cases 

XChICP 9 LOCATIONS Op CHORONISE CONTROL POINTS 
XETACP r LOCATIONS OF SPaNWISE CONTOROL POINTS 

TF( .NOT, CHECK) GO TO 60 
IFdOt.EQ.lDDn GO TO 46 

iF(eATCM) CALL ST0P2(UC0d rC FILE ANO AIM FILE OO NOT 

jHAVE Matching lo • iFloatciod) 

Pause • bc file and aim file do not have matching id • 
48 IF(I005,EU,I05) GO TO 60 
IF(BaTCH)CaLL ST0P?(UC0, 

1» 1DD3 ODES NOT Match bc file value SPL0AT(ID5)) 

Pause »idoi does not hatch bc file value» 

check THaT CHORDHISE aNO SPaNwISE LOCATIONS MaTCH ON 
THO FILES 
60 DO 61 1 *1 ,PP 
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IF(CHlC»*(I)tEO,XCHICPfI)) GO TO 61 
IF(BATCH) caul 5T0P2 

l(U6f» CHOROWISE CONTROL POINTS ON FILES 00 NOT MATCH 
aFLOATfl)) 

Pause ^chordwise control points do not match^ 
fc1 continue 

DO 62 I » 1 ,NMP 

lP(€TACPfI),EQ,XCTACP(n) GO TO 62 

IF(BATCH) CALL ST0P2(U6# '3PANNISE CONTROL POINTS DO NOT 
1 Match on ? files aT index »,ELOAT(n) 

Pause tsPANwisg control points do not match* 

62 continue 

initialhe variables 
KJUMP • 2 FOR symmetric WlNG 

B 1 for UNSYmmETRIc ring 

IFfSYM) KJUMP • ? 

IF(ASYM) KJUMP B I 
KK? B KK 
NMAX? B NmAx 
nROWS « PP 
U B 0 
li • 0 
1 5 B 0 
PPnEN b PP 
nTYPE b 0 
HMNEN s MM 

middle b M0D(MM#2)bMmP 
IF (SYM) GO TO 6S 

IF (Modes. gt.O) go to to 
if(8atch) Call sroP2fU6, 

X I UNSYH MING BUT NOT ALL MODES ARC AVAILABLE >i 
2 FLOAT(MOOe5)) 

Pause «unsym ming but not all modes are available « 

6B IF(CONV) WR|TE(UrO#BaO) 

read in value for NTYPE 
0 8 symmetric cases 
-1 B antisymmetric cases 
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C 1 s BOTH SYMH&TRIC AhD ANTISYMMETRIC CASES 

C fOR an UNSYMMtTRIC WING NTYPE « 0 

REAOfUCl f BIO) NTYPE 
IFrnOOES^GT.O) CO TO 85 
IF (MOOES, EQ,0 , and, NTyPE,EO, 0) GO TO 75 
IF (MOOES, LT,0, AND, NTYPE,LT,0) GO TO flO 
IP (MOOES, EQ.O, and, nTYPE.GT, 0) GO Tt) b8 
IF (MOOES, LT,0, and, NTYPE, GT.O) GO TO 69 
IF(BaTCH, AND, NTYPE, LT, 0) CALL STOP^cUb, 

II aNTi-SYMMETRjC MODES ARE NOT AVAILABLE W 
^FLOaT(MOOES)) 

IFfBATCH. and, NTYPE, GE, 0) CALL 
iStOPaCUb,! SYMMETRIC MOOES aPE NOT aVaILaBLE 
2FL0AT(M00Esn 
IF(NTYPE.LT,0) 

JPaUSE f aNTI-SYMMETRIC mooes are not available* 
if(ntyp£,ge,o) Pause 

pSYMMETRIC MOOES a«E NOT AVAILABLE* 

GO TO 65 

68 NTYPE a 0 
wRITE(UCOj88n 
GO To 75 

69 NTYPE e »1 
WRITE<UC0#862) 

GO to 60 

INOOEL 1$ AN array CONTAINING INDICES OP COLUMNS OR 

rops to re deleted 

indrow is an array initially filled WITH indices of all 

ROWS FOR THE PARTICULAR TYPE OF WING, 
iNocoL IS an array initially filled with indices 
OF all columns for The Particular case’, 
after Each deletion a call is m^de to subroutine krunch 
after which InOROw contains indices of rows to bE 
retained or indcol contains indices of columns to 
BE retained, 

NRECD s number of LOGICAL RECORDS 
NCOLSB s number of columns before deletion of 
unwanted columns 



ririn or»o no oooonnor* 


NCOLSA ■ NUMetR OF columns after deletion of 
unwanted columns, 

NROWSe B number of rows before deletion of unwanted 
Rons 

NROWSA s number of RQWSA AFTER DELETION OF UNWANTED 
Rows, 

UNSVMMETRIf WING 
70 NXR B PF * MM 
NCOLSBbNMAX * KK 
NRECDS B MM 
GO TO 120 

symmetric boundary conditions only 

7% NXR 8 PP ♦ MMP 

NcOLSB 8 NMAX *(fKK^n/2 ) 

NRECDS 8 MMP 
NXC 8 NCOLSB 
GO TO 


SYMMETRIC WING WITH aNTUSYMMETrIC B, C, ONly 
SO NXH 8 PP * MMPA 
NRECDS 8 HmPA 
NXC B NMAX * fKK/2) 

NCOLSB B NxC 
GO TO 121 

symmetric Cases for symmetric wing with both 
symmetric aNO aNTUSYMMETRIC boundary conditions 
IF(NTVPF,LT,0> go to 
NXC 8 NMAX 8((KK4n/?) 

NCOLSB 8 NMAX t KK 
NXR g PP * MMP 
NRECDS 8 MMP 
lY 8 0 

DO 90 I 8 1 ,KK,^ 

L « f l*n * NMAX ♦ I 
K » L ♦ NMAX » 1 
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DO 90 J « L,K 
IX c IX 4 1 
90 I^JOCOL(IX) » J 
r,0 TO l2S 

*NTI.SYMMaT«lC CaSIS FOR A SYWMCTRIC WING WITH j^OTH 
symmetric and antusymmetric boundary conditions 

95 NXC « NMAX * f KK /B) 

NROWSB « PP * MMP 

space past The symmetric boundary conditions as they 
are not needed 

00 9h I ■ l.NSYM 

96 RtADfUS) 

NCOUSB c NMAX * KK 
IX a 0 

00 97 I c 2 ,KK,p 

1 a fl*n * NMAX ♦ I 
K a I ^ NMAX • 1 
DO 97 J s L,K 
IX a IX 4 1 

97 INDCOLtlx) • J 
NXR a MMPA « PP 
NRECDS a MmPA 
GO To 

120 NXC a NMAX 4 KK 

121 00 122 I » 1,NXC 

122 INoCOLfI) a I 

123 00 12« I a l.NXR 
129 INDRON(I) a X 

NROWSe a NXR 
N2 a NxC 
N5 a NxR 
NCOtSA a NXC 
nROwSA a NXR 

C.*.,t^SYM2 ANO NASYM2 ARE THE NUMrERS OF SYMMETRIC AND ANTU 


SYMMETRIC SOLUTIONS TO BE OBTaInED aND BCS? 

C AMO BCaS? DENOTED THE TYPES OF CASES, 

C 

NSYM2«NSym 
NASYM^ aNASYM 
DO ISO I«I»|0 
BCS2(I)a0C5(D 
ISO 6CA$2(naBCAS(n 

IF (UN5YM,Nr,0 ,OR. nTYPf,CT,0) GO TO 150 

if(ntype,lt,0) go to l<*0 

NASyhZsO 
DO 132 lal.lO 
BCAS2(l)a. false, 

132 CONTINUE 
GO TO l5o 
lao N$YH2ao 

DC loa laiftO 
BCS2(I)a,EALSE, 

„ ia2 CONTINUE 
150 CONTINUE 
C 

NCOm 8 0 

IFfCOM,EQ, KNE*) GO TO 200 
C 

c initialize commands 

153 STORE a .TRUE, 

COMSUM a ,FAlSF. 

POINTS « .false, 

RH5 s .false. 

PSOL a .FALSE, 

CONTIN s CONV 
PIND 8 .false. 

200 CONTINUE 

NCOM a NCOm ♦ 1 

IF (NcOM.GT.NCOMK) call ST0P2rU6, 

1» TOO many commands float (NCOm) ) 

C 

IF (cONV) MHlTE(UeO»0«5) 

REAO(UC1#BoS) (CAR0IM(I,NC0M) , ,a) 


COM • CAROIM(UNCOM) 
aiO CONTINUE 


IF 

(COM 

, EQ • 

Hals) 

GO 

TO 

225 

IF 

(COM 


NEWS) 

GO 

TO 

It 

IF 

(COM 

• to. 

knes) 

CO 

TO 

17 

IF 

(COM 

1 E Q . 

OCMf) 

GO 

TO 

sao 

IF 

(COM 

• to. 

DSM$) 

GO 

to 

360 

IF 

(COM 

• to, 

OCRs) 

GO 

TO 

3T0 

IF 

(COM 

• to. 

DSPS) 

60 

TO 

3Q0 

IF 

(COM 

• to. 

PPPS) 

GO 

TO 

396 

IF 

(COM 

,EQ, 

CONS) 

GO 

TO 

235 

IF 

(COM 

• to. 

stas) 

GO 

TO 

aoo 

IF 

(COM 

.to. 

NSTS) 

GO 

TO 

211 

IF 

(rOM 

.to, 

STOs) 

GO 

TO 

212 

IF 

(COM 

• to. 

SUMS) 

GO 

TO 

213 

IF 

(COM 

• to. 

NSUs) 

GO 

TO 

2ia 

IF 

(COM 

,to. 

PCPS) 

GO 

TO 

215 

IF 

(COM 

• to. 

NPCI) 

GO 

TO 

216 

IF 

(COH 

.tQ, 

RHSs) 

Go 

To 

217 

IF 

(COM 

• to. 

NRHS) 

GO 

TO 

216 

IF 

(COM 

• to, 

PSOl) 

GO 

TO 

2l<9 

IF 

(f OM 

.to. 

NPs$) 

GO 

Td 

220 

IF 

(COM 

.to. 

PINs) 

GO 

TO 

221 

IF 

(COM 

• to. 

NPls) 

GO 

TO 

222 

IF 

(COm 

• to. 

SETs) 

GO 

To 

153 

IF 

(COM 

• to. 

TSSs) 

GO 

TO 

223 

IF 

(COM 

• to. 

RPMI) 

GO 

TO 

226 

IF 

(COM 

• to. 

RrPs) 

GO 

TO 

232 


WRITE fUCOf 650) (CaROIm(I#NCOm) I l»l |i*) 
If (, not,cOntin) stop 

NCOH ■ NCOm m\ 

GO TO 200 

21i STORE s .false, 
loaio 
GO TO 200 

212 STORE 8 .true. 

GO TO 200 

213 COMSUM e .TRUE, 
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no TO 200 

i\U COMSUM 8 ,FAUSI:, 

GO TO 200 

2JS POINTS 8 ,TKUf:, 

GO TO 200 

216 POINTS a .false, 

GO TO 200 

217 RMS 8 ,TpUF, 

GO TO 200 

21ft RMS • .FALSE, 

GO TO 200 

210 PSOL # .TRUE, 

60 TO 200 

220 PSOL M .FALSE, 

GO TO 200 

221 PINO 8 .true. 

GO TO 200 

222 find 8 .false, 

^ GO TO 200 

o» C 

223 CONTINUE 

Tss command is for The ames* version ONLY. IT allows commands 
TO rE given to the operating system from the program, 

REaD(UCI,606) TSSCOM 
call OBEY(fto,TSSCOM) 

GO TO 200 


Halt 

225 IF(OPENb) ENnFILE U6 

IF(0PEN6) Call ObEYMSiISHRELEaSE FT0<*F00l ) 
WPITE(UCOr855) 

GO TO 999 
C 

226 CONTINUE 
C 

c rpm-.-restore pressure modes 
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KK^ s KK 
NMAX? c HHkH 

IF (HOOES,GT*0 ,AND, SYM) GO TO 226 
DO 227 I®l,NCOLSe 
227 lNDCOLcI)»I 

NCOtSAsNCOlSB 

N2«NC0l5B 

r,o TO aoo 

226 CONTINUE 
C 

C MO«E NODES A^t AVAILABLE THaN REQUIRED, I,E,, SYMMETRIC MING NITH roTH 

C.,,, .SYMMETRIC aNO ANTISYMMETRIC MODES AVAILABLE' 

C 

KL»I 

IF fNTYPf.tT.O) KL*2 

NCOLSA at 0 

DO 22R K«KLfKK»2 

I«PP*(K*n 

DO 22R P«l,RP 

NCOtSAcNCOLSA^i 

INDcOL(NCOlSA) 8 P+! 

229 continue 
N28NC0USA 
GO TO 200 
C 

2S2 CONTINUE 
C 

C RCP— RESTORE all CONTROL POINTS 

C 

Ll »0 

L2»0 

lJ*0 

LIST! (1) * 0 
L1ST2(1) s 0 
PPNfNsPP 
MHNEusMM 
NROwSAsNROwSB 


C 
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antisymmetric Cases ON^Y the control ROINTS aT the CENTERlINE 
,,..,HaVE to be OELETEO; control points wiu BE AT THE CENTERLINE 
HM is 000 (H0D(MM,2)aj), 

lEf M0 DES.lt «0) NR0w5ABNR0WSA»M0r>(MM,2)*PP 
N3 ■ NROhSA 
DO 233 Ut .NROi^SA 
33 lN0R0Wrn«I 
60 TO 200 

CONTINUE 

23S CONTIN 8 ,TRUE. 

60 TO 200 

OcM — —•delete cHOROi^ISE MOpES 
NMa)(2 is the new NlIMpER OF CHOROWISE MODES 
340 IF(CONV) WRlTEfUCO.fiSS) NMAX;»,NmAX 
READ(UCl«aiO) NMAX2 
NHAX2 s HINOfNHAXf NHa)^2) 

IF(nmax2.eq.nmax^ go to ?00 

tF{NMAX2,GT.O) go to 3RR 

IF(BaTCN) call 5T0P2(U6n INCORRECT VaLUE FOR NMAX2 »# 
1FL0AT(NMaX2) ) 

WPITE(UcO>867) 

NCOH 8 NCOh • 1 
GO TO 200 
344 J « 0 

lADD 8 0 

IHAX 8 NMAX2 ♦ 1 
346 DO 346 1 a IMAX.NHAX 
J s J* I 
L a I f lAOO 
INODELfJ) 8 L 

IF(IN0D£L(J),£Q.NC0LSH) GO TO 3S0 
348 CONTINUE 

UDD 8 lADD t NHAX 
C;0 TO 3^6 

350 call K«unCH( J, IN0DEL,N2, INOCoL#NCoL8A) 

Np 8 NCOlSA 
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GO TO 200 


OSM---DElETt SFaNWISE MOOES 
KK^ IS The new NUMbEF of SP 4 NWISE MOOES, all 
columns to The fight op kk? * nma)( mill be deleted, 

360 IF(CONV) NMlTEfUCOfSTO) KK2^KK 
REAOcUCIfSjo) 

KK2 s HiN0(KK|KK2) 

IF (,KK2,fO.KK) go to 200 
IF{KK2.GT,0) GO TO 364 
IF(CONV) GO TO 36? 

Call SToP2CU6#« INCoRBCCT value FqF i»FLoAT(KK2)) 

36? MRITE(UCO#e7o) 

NCOM t NCUm • 1 
GO to 200 

364 IF Rhodes, GT. o) ix m kk2 * nhax ♦ i 

IF(MODE$,EO,0) It s ((KK?fn/2)*NHAX4i 
IF(MODES,LT,0) LX • (KK?/?)*NMAX^i 

J 9 0 

00 367 I c Lx,NC0LSB 
J 8 J 4 1 
INODEL(tJ) m I 
367 CONTINUE 
GO TO 35o 

DcF-«*OELETE CHOROWISE POINTS 
POEL IS aN INTEGER aRRaY OF NOT MORE THaN 
CHOROWISE POINTS TO «E DELETED. MUST qE 

IN ascending order; 

370 IF(CONV) WRITerUcO,87S) 

REaOcUCIi 076) (PDEL(l)»l8l#l7) 

IFCPOELU) lEOfO) GO TO 200 
DO 373 I « 1.17 
IF(PdEL(T),EO,0> go TO 372 
Lt 5 LI 4 1 
LISTKLU 8 PdEL(I) 

C LI 8 T 1 IS AN array containing THE INDICES OF DELETED 
C CHOROWISE POINTS, 

60 TO 373 


37? LAST B I - I 
GO TO i57« 

373 CONTINue 
I A5T»17 

37<i PPN|W a PHNEW . t.AST 
Nl • 0 

DO 376 I » l^NktCDS 
DO 376 K a I, LAST 
Nt « ♦ J 

IPOIMT ■ PP * fl-n ♦ POEtfK) 

376 IN006L(^»l)» IPOINT 

38a CALL KRUNCMCNj # lNDDiL»N3, IN0R0W,N«DN5A) 

N3 • NKOWSA 
r.O TO 200 
C 

C DSP— delete SPaNWISE points 
c MOEL IS aN integer array OF NOT MORE THAN |7 
C SPaNWISE POINTS TO be DELETED IN ASCENDING ORDER, 
_ 390 iFfCONV) WRlTEfUCD,877) 

o C TO delete these POINTS INCREMENTS OF PP ARE 

C added TO reach point (NRECD.n TIMES AND THE 

C results stored In THe INdDEL aRRaY, SUBROUTINE 

c KRUNCH IS Then called to delete These points. 
REaD(UC 1#876) (MDEUniI»l#l7) 

IF(MDELn) ,£0,0) GO TO ?00 
N I a 0 

DO 394 I « 1.17 
IH » MOEL(I) 

IF(IX,fO.O) go to 384 
1F(3YM*AND.xX,GT,MjODLE)GO TO 594 
HMNEW a MHNEW • 1 

1F(IX,NE, middle, aNO.SYM) HMNER a MMNEW • t 
L? * L? ♦ I 
LlST?a?) « 

LX a (IX. n * PP ♦ 1 
JX • IX ♦ PP 
00 39? J ■ LX,JX 
Nt a N! 4 I 
39? INDDEL(Nl) a j 


on o 


COMTlNijE 
(50 TO nea 

opp— otitTE pa«t:cula« point 
read indices of points to bE deleted 

39B IFfCONV) HRITEcUC0#876) 

PEAO(UCX,StO) P^M 

LS • IJ ♦ 1 
LISTinrLJ) « P 
UST3(2,L3) « M 
iNDOFin ) * fMi^n * pp f p 
NJ • I 
CO TO X8tt 

c 

c start 

UOO IF (NROwSa.GE.NCOLSa) go to aio 

IFfBATCH) CALL ST0P2(U6^ HORp UNKNOWNS THaN fOUaTIONS 

IFLOaT(NROWSa) ) 

Pause «more unknowns than equaTionsi 
GO TO 200 

c 

R20 iFf, NOT, STORE) GO TO aSO 
C 

CALL OBEY f22f2RH0DEF FT09 FOO W , IdF RE ) 

C 

rewind U20 

REaD{U2o) lOAf I0B»10C«IDa 
toa ■ tD4 ^ I 

rewind uao 

WRlTE(U2fl)lDAf IDBf IDC# 10a 

endfile u2o 

call 0BEY( 16, 16HRELEASE FT09F001) 

IF (OOISK.NE.O) WRITE(UC0#932) ID« 

C 

C WRITE heading information 

430 write (U6, 866) ( TI TLE ( I ) # I« 1 1 NT ITL ) , lOl • 102 , 103, lOa , 

I JJ,UNSYM,M0dES,NTYPE,PPNfW,MMNf;w,NMaX 2,KK2^ 
2DELTAO,MACH,<BCS(I),I*l,lO),(BCASn),I*l,lO) 

WRITE(U6,8i2) 


WRITE summary OE COMMaNOS 
IP (t^^0T,cClM5UM) GO TO aao 
WRiTECUbrfiSn 
DO «^5 J « t,NCOM 
WRITE (U6f807) ( C ARO I M ( I , J ) , 1 a l ^ ij ) 
ai5 CONTINUE 

SORT OEIETEO indices 
« a0 iFf , NOT, POINTS) GO TO «77 
lEaifLEfl) GO TO aaa 
IX » Ut 
K 8 1 1 • 1 
DO aaa j k i ,k 
IX f Tx • 1 
DO L « 1 , IX 

I 8 L ♦ 1 

IFaiSTJ(U),LE,USTlfI)) CO TO ua? 
HOLD « LISTi ft) 

ItSTl (1) 8 LISTI f n 
LISTun 8 HOLD 
RR? CONTINUE 

««4 IE a2,LE,n GO TO «7o 
IX a L? 
k z Li m I 
DO 446 J a 1,K 
IX 8 IX . 1 
DO 446 L a | , I x 
I a L ♦ 1 

lF(LIST2(L),LEai5T2f in GO TO 446 
HOLD a LlST^fL) 

LiST?(L) • LiSTafi) 

LIST?(I) a MOLD 
446 CONTINUE 

WRITE OUT fHOROWISE CONTROL POINTS 
470 WRlT?(U6r9lO) 

IX a 1 


n r» o o o r» o r* 


lF(UlSTUIi(),NF.n GO TO 472 
IX s tX * 1 
GO TO 475 

472 WRlTefU6,906) CHICP(I) 

47'^ CONTINUE 

C WRITE OUT SPaNWJSE control POINTS 
WRITE(U6,9i2) 

IX s i 

DO 475 I * l.NRECOS 
1F(LI5T2(IX) .NE,I) GO TO 474 
IX ■ IX f t 
GO TO 475 

474 WRlTE(U6,9o6) ETACP(I) 

475 continue 

ir(LS,ro,o) GO to 477 

c 

c write out DtLETEO CONTROL POINTS 

WRITE(U6,9|4) 

DO 476 I » 1,LJ 
J * LISTSOfl) 

K s LISn(?,l) 

476 WRITe(U6,915) CHICP ( J ) # f T ACP ( K ) 

C 

C read »OUNDaRY condition file 

477 HERO « 0 
IF(SKIP) GO TO 486 

lENO s NUMBER QP RIGHT HAND SIDES 

IF(NTYPE,GE,0) IENO b NSYM 
IFfNTYPE,LT,0) IENO « NASYm 
jFfASYM) jEND * NSYH ♦ NASyH 
ISTOP s UNO 

lFfNTYPE,GT,0) iSTOP b NSYM ♦ NaSVM 
both sets OF boundary CONDITIONS A«E READ IN aT THE 

Same time if both types of solutions are required, 

This eliminates ThE need to rewind ThE boundary condition 
file for subsequent Cases op the same configuration. 


IX « NHOwSB 
00 a79 1 » 1,IST0P 
IPfNTYP^.NE, n GO TO <i79 
1F(1,GT,WSYM) IX 8 MMPA * PP 
U79 READ(U8)(ALFA(K^I),K«1,1X) 
a«0 IFC.NOT.RHS) GO TO «B6 

WRITE boundary CONDITONS 
WRITF(Ue,# 8 <? 8 ) 

00 aea n « uiend 

1 8 IZfcRO ♦ II 
WPITECU6,697) I 

J * I 
K 8 PP 

DO asi M s I ,NR£CDS 
WRITE(U6,899) (ALFA(U, n,L»J,K) 

J » J ♦ PP 
46*^ Kc K ^ PP 

WRITE(Ue», 8 l 2 ^ 

«6a continue 

UB6 IFf ,NO t,PXNO) go to 

WRITE OUT INOcOL. INpHOW 

WRITf(U6,890) 

WRlTf (U6,692) (INDCOL(J) ^UsUNCOLSA) 

WRITf (U 6,d9a) 

WRITF {U6,89;>) (iNOROWf J) » J»i #NROWSA) 

CALL equation solving SUrROUTINE 
09) NkTP * NCOlSA 4 lEND 

write fU0#933) NROwSAf NCO lSA 

Call solvit (NC0LSB#NC01 Sa^NROWSB, NROWSa, iend^reaoa, 

I IFHT,Ui 1 ,INOcOL, INDROW,NR£cDS, InVERS, A t ALF a, IZERO,Nh#NX 
aNAlNV, ATA, AR# IL#X, AlNV,RH$REStOETt lEPPf NKTPtNROWS) 

IF f IERR,NE,0) CALL STOP^rUbt 
p ERROR IN tnUAtjON SOLVERt lERR IS f, 

2FL0AT(IfRR)) 


n r» o n o r» n 


WRne(U6,695) d£t 


SOLVIT OOtS NOT calculate HESIOUaLS CYET), 

WRITE CUCO#90^) 

DO «95 I s I • IEND 
WW1TE(UCO#«>SO) If RMSRfSfl) 
i|Q5 CONTINUE 

lrf,NOT,PSOL) GO TO S12 
C 

C PRINT SOLUTIONS FIVE CASES AT A TIME 

IFrIND,NE,n CO TO 5o3 
WRjTE(U6f902) 

GO TO 

5o5 IF(ASYM)wRITE(U6,9oo) 

If (SYM,AND,NTYPEtCE*0>‘^^I^E(US#90<n 
IF($YM,ANO,NTYPE,LT,O)W«lTE(U6f902) 

C 

Soa CONTINUE 

DO Si 1 JJJ» 1 f IEN0,5 
J9 B MlNPCJJj4«f l^NDl 
WRlTE(USf9on (CASE, JKf JK.JJJ, ja^ 

I C OUNT B 0 

DO 5o5 II B l^NCOLSA 
IFClCOUNT.rOf 0) WRlTprU6,8l2) 

ICOUNT m ICOUNT ♦ I 
IFflCOUNT E0 ,NMaX 2> ICOUNT a o 
SOS write (US#90T) (Ilf (itdif JK) f JKaJJJf ja) ) 

511 CONTINUE 

512 1F(, NOT, STORf) GO TO 5tS 
C 

C STORE results 

c IF both symmetric and antiysmmetric solutions are to 

C STORED BYPaSS OPENING FILE AND WRITING TITLE INFORMATION, 

iFtlND.EO.n gO to Sl3 

C.*..t8NKFIL IS FOR aMES» TSS VERSION, rNKFIL ISSUES CONTROL CaRoS TO The TSS 
C, , .operating system to define (ODEF) the SOLUTION FILE, 

C 



ooooon r* o o r> r> n r> r» r» 


CtLL SNKFIL(I0)>10?^ID3>t0«) 

WRITE(U12)IDl.I02,103,IO<t,TirLE,SYM,N5YM2,N*8YM2, 
uecs?(n, 1*1, 10) , {BCAsjd), HI . JO) ,nm*x?,kk2,jj, 
2PPNEw,MMNEw,n«ONSa,OElT*0,EPS,m*CM 
5lJ DO 5l« I » I.IEND 
Si a mrITe (U 12) (X ( J, ! ), J*l fNCOLSA) 

5l« IF(NTYPE,LF,0) GO TO 550 

CASE WHERE THE WING 18 SYMMETRIC aND SOtUTIONS A«E 
TO be obtained for both symmetric and aNTI-SYMMETRIC 

BOUNDARY CONDITIONS 
IFdND.EQ.l) GO TO 5«0 

add nmax to Each element of fiNDCOU to obtain 
antisymmetric mooes 

DO 522 I « I.NCOISA 
522 lNOCOt.fI> t INOCOltI) ♦ NmAX 

lENO ■ NaSyM 
IND c { 

IZERO ■ NSym 

IF(MOO(KK2,2),NE,0] NCOLSA 1 NCOISA«NHaXj 
nCOlSB remains the Same 
IF (MMP.EO.MMRa) go to 526 

deleting control points on centerline, 

A DIRECT count of POINTS ON THE cENTERLlNE IS 
NECESSaRY since the DPPOINT command may have been 
USED, NUMCL ■ number currently ON THE CENTERLINE, 
NUHCl ■ ft 

JMlN « NROWSB«PP*l 
IMIN * MAXo(NROWSA»PP+} , I 1 
DO 525 J * JMjN,NROWSB 
DO 523 1 ■ 1M1N,NR0WSA 
lF(J,£a,INpROW(I)) 00 TO 5?fl 
523 continue 
GO TO 525 

5ju NUMCL » NUMCL ♦ 1 



ono n r> r» o o o r» ooo r>r*o 


525 IHIN « MINOtlMlN^J ,NR0WS4) 
nKOWSA 3 NRO^SA « NU^Cl 
nROMSB « nROMSB • PR 
NRECDS b NRECDS • 1 

52«» continue 

backspace the Influence matrix file to the beginning 
OF the influence Matrix, 
remind un 
read (uin 

IF skip b .true, bypass READING OF BOUnDaRT 
CONDITION FILE, 

IF(SKIP) GO TO <lBS 
GO TO aBfl 
540 CONTINUE 

NON that the antisymmetric solutions Have been 

obtained AND (IF REOUlREO) STORED# THE 

affecteD Variables kiu- be reset to valufs appropriaTf 

FOR further solutions TO THE SYMMETRIC CASES 

lENO b NSYm 

I NO • 0 

1P(MOD(KK2#2).NE.O) NCOlSA • NCOlSA ♦ NMAX2 
DO 542 I B 1, NCOLSA 
542 INDCOL(l) » INDCOL(I) - NMax 
IF(MMP,f«,MMPa) go to 546 

NROnSA b NROnSA 4 NUMCL 

NRONSB b NRONSB 4 PP 

NRECoS B NReCdS ♦ I 

546 CONTINUE 
550 CONTINUE 


IF(STORE) ENoFILE U12 
IF(STO«E) WRlTFrU6#685) 

aT This point the soLUTioNrS) haVE been obtained 

AND(IF REOUIREO) STORED ON THE SOLUTION FILF, 



o r» r> r> r» r» 


IF MORE SOLUTIONS aRE TO BE OBTAINED WITH THE CURRENT 

influence Matrix, uij rill have to be rewound and 
SPACED TO the beginning OF THf INFLUENCE MATRIX, 

THE PROGRAM ASSUMES THIS IS THE CaSE UNLESS THE NEXT 
command is Halt, new, or knew 

WRITE(U6f 868) 

IF(CONV) Wr1TE(UC0^8«5) 

NCOM « 1 

READ(UCI# 80 S)(CAR 0 IMfI,NC 0 M),lal,a^ 

COM a CARDlMf I ,NCOM) 


IF 

(COM 

4 E U t 

haLs) 

GO 

TO 

225 

IF 

(COM 

• eq, 

NEWS) 

GO 

TO 

17 

IF 

(COM 

• to, 

KNEl) 

GO 

TO 

(7 


SKIP ■ ,TRUt. 

NCOm • t 

r 

rewind Ull 
^ read (Ult) 

00 write (U6^8t2) 

GO TO 2l0 
C 

7R9 FORMATfUn 

8ol FORMaTCJH , CENTER BaTCH»/) 

805 fORMaT(RaS) 

808 FORMATr^OAa) 

8o7 FORMaThH ,«a5) 

BtO f0RmaT(16IS) 

812 FORMaTHH ) 

8iS FORMaTCjH ,»ENTER OOISKt) 

820 fORMaT(1h #*eNTeR lOOl , 1002# ID05 • ) 

aao formaToh #iENTEr value for ntype# o por symmetric# »#ix# 
!»•! FOR aNT1»SYMMETR1C# »# ’ I FOR h0TH» ) 

805 FORmaTIIH , !♦» ) 

850 FORMaTOH # tUNRECOGNIZED COMMAND f,ttAS ) 

65S FORMaTUH #<PROGRaM TERMINATED BY HaLT COMMaNQ'/) 

865 FORmaTCIH #iEnTER nEW nUmHER OF CHORDWISE MODESt# 
llX,iOLO VALUE ISi#I5, • ORIGINAL VALUE WASi#I2 ) 

867 FORMaTIjH #»tNCORRECT VALUE FOR NmaX?, aT ♦ SiGNt, 


L .. 



11X,»REENTER COMMaND OCM t ) 

070 FORmATUH NEw mAXImUM NUhBER of SRANi^ISE MODES', 

nx,tOUD maximum is *, 12/* original, value Was 1,12/) 

87a formatcih .umcorrect Value for kk?, aT ♦ sign t, 
pRHENTER command OSm* ) 

875 FORMaTUH ,»ENTER cMORDWiSE POINTS TO rE DELETED *, 

iix,*i7 maximum in ascending order* ) 

878 FORMATnnS) 

877 FORMaTUN ,*ENTER SPaNWISE CONTROL POINTS TO hE * , 
llX, iOELETED, \7 maximum ASCENDING ORDER* ) 

878 FORMaKIH , ‘enter INDICES Of THf PARTICULAR *,IX, 

PPOINTS TO 0E DELETE0—P,M* ) 

881 F0RMATC//»0ANTISYMMETHIC MOOES a«E NOT AVAILABLE.*, 

12X, iPROCEEDJNG with JUST SYMKeTrIc MODfS,*//) 

882 FORMaTC/Z'OSYMMETRIC MOOES ARE NOT AV A I L ABLE , » , 2X , 
1*PR0cEE0ING WITH JUST ANTISYMMETRIC MOOES*//) 

085 FORMaT(IhO, ‘SOLUTION FOUND AND STORED*,/) 

686 FORMaTUHI/* determination OF BNK COEFFICIENTS*/ 

^ j ♦ * / 

^ nY,20A^i//* IDl » ',13/' In^ » 

1* IDl • *,I3/' loa s *fIl/»* JJ • 

UNSVM s *#I3/,* MOOES i *|l3,/» nTYPE x », 

5 lj/* PPNEW ■ *, 17 /* MMNEW X *, 15 /’ NMAY? » *, 

All/* KK2 m *,il/* OELTAO = *,F6,a/* mACh « *f 

5F8,a/* BCS X »,10U/* BCAS x *,10l1/) 

887 FORMaTCIH , iSUMMaRY OF COMMANDS*,/) 

888 FORMaTC//' finished^/) 

8R0 FORMaTCIHO, * INDICES OF COLUMNS TO BE RETAINED*/) 

89? FORHaT(7(iX,!3,1X)) 

B9a format (IHO, * IND lCfS OF ROMS TQ BE ReTaINeD*,/) 

895 FORMAT!/* THE DETERMINANT OF ( A-TRANSPOSE) *A r •, 

IEIS.8/) 

897 FORMaTCIHO, *CaSE *,I2/) 

898 FOrMATUHO, 'BOUNOArV CONDITIONS AT OrIGINaL CONTROL POINTS*, 
I/* (L,F, TO T,E.. RIGHT TIP TO LEFT TIPi*/) 

099 F0RmaT(IH0,8(F10.5,1)O) 

900 FOhMaTMHO, 'SOLUTION FOR UNSYMMETrIC CaSES*/) 

901 format (IhO,2X,5(5X,A«, II, IX)) 

902 FORMaTCIHO, 'SOLUTION FOR aNT I*SyMME TR I C CASES*/) 



C FORMaK/* CaSE» rSX, tRH8 RESIOUAU'/) 

90^ Fo»HaT(JMO, fSOLUTION FOR SYMMETRIC CASES’/) 

906 format flH ,M0,5) 

907 FORMATftH lX,SfM0.7,3X)3 

910 rORMaT (iHOr »CHOROHlS|r cONTROt POINTS UStOV) 

9|2 FORMATUHOr 'SP aNNISE CONTROL POINTS USED* ,/) 

9ia FO»HaT(iHo, ’SPECIFIC POINTS DELETED’/’ CHlcP’^ 
llOX^ *ETACP'/) 

9t5 FORMATriH r FB ,S # 9x , P 8 p S ) 

9J0 F0RHATf^X#I^p7X*FJ0,6) 

931 FORMaTC ' OUTPUT jS ON ,,, OUTPUT , SOLVR , N » p j I , 3 h t ,. ) 

932 FOHMATP loa«’,J5) 

933 FORMaT(» number OF CONTROL POINTS *WIS/ 

1 » number of pressure MOOfS «»#I5/) 

999 STOP 
End 

subroutine 5QLVlT(NC0L5BpNC0LSAf NRONS8pNROMSA,NRHSS^READAp IFMT, 
iNUNITpINOcOL, INOROWpNRtCD, INVER8»AiBp IZEROpNB^NX^NAlNV, 

» 2ATA,ARp IL#X,AlNVpRMSR£S,DET,lERRfNKTP,NROwS) 

® C 

c 

c THIS SUBROUTINE SOLVES a SET OF SIMULTANEOUS EQUATIONS 
C (A)(X)e(B) AND/OR determines THE INVERSE MaTRIX FOR A SQUARE OR 
C rectangular (A) matrix IN THE LEAST SQUaHeS SENSe USING 
C HOUSEHOLOERiS method fUNlTARY TRl aNGULARI^aTIOn OF a N0N« 

C SYMMETRIC MATRlXf U, ASSOC, COMP, MACH, 1958). 
c THE (A) matrix Can either be Input through the argument [ ist 
c or read from aR external device, in either Case the program has 
c the capability of striking out rows aNO/OR columns before deter. 
c MINING The solution ANO/OR The inverse, 

C 

c INPUT’ 

c — — 

c NCOlSB m number Op COlS, REpORE DELETION Op UNWANTED COLS'. 

C NC0L5A m NUMBER OF COlS. AFTER OELFTION OF UNWANTED COLS. 

C NROW8B » NUMBER OF ROWS BEFORE DELETION OF UNWANTED ROWS. 

C nROwSA s number of ROWS AFTER DELETION OF UNWANTED ROWS, 

C NRHSS s number OF RIGHT HAND SIDES, 

C REAOA 9 LOGICAL VARIABLE WHOSE TRUTH CAUSES THE (A) MaTRTX 


r»r»o or>fir»r»r»oonr»oonor»rir» c*»r» r»r>r»r>nr»r»ori rir»ooorir»o' 


rFMT 


NUNH s 

( iNOCOUc 
(INDMOw), 
nROnS c 


NKTP 


(6) 

' a 

Nr 

a 

NX 

a 

NAXNV 

a 

NRECD 

a 

HERO 



TO BE read from an external device, 

FORMAT indicator (RE A0A«, TRUE , ) . 

0 means unformatted 

1 MEANS format raFlO,0) 

? means format (SFIS.O) 

% means free format («F) (IF supported by COMPILER), 
anything else causes a return MITH IERR«r. 

FORTRAN UNIT NUMBER FROM WHICH TO READ (A) (REAOA* 

•True,), 

indices of columns TO 8E RETajNEO, 

indices of rows to be retained, 

the number of rows stored in Each logical record 

OR THE number of ROWS YOU WISH TO PROCESS aT ONE TIME 
(READA«,TRUE,) or the first DIMENSION OF (A) IN THE 
program WHICH DETERMINED (A) («EADA»,f ALSE , ) , 

THE number of COLUMNS OF aTa. SHOULD 0E SET IN CaLLING 

program to NCLOSa * NRHSS FOR I NVERSa , F ALSE , 0« TO Two 

times N^OLSA ♦ NRHSS FOR I NVcRSa , TRUg , 

THE SET OF RIGHT HaNO SIDES BEFORE STRIKING OUT UN* 

wanted rows (b) will not bE destroyed, 

THE FIRST dimension OF (b) IN THE PROGRAM wHIcH DETER* 

MINED fB), 

THE FIRST dimension OF fX) IN THE PROGRAM IN WHICH 
(X) WILL BE USED, 

ThE FIRST Dimension of (ainv) in the prooRam in which 
(ainv) will bE used, 

number of logical records or number of times YOU WISH TO GO 
THROUGH REaO loop, NRECD TIMES NROWS MUST BE EQUAL TO 

nrowsb Plus o to rows-i 

COLUMN OF B ARRAY WHERE RIGHT HAND SIDES FOR 

THIS Case begin 


OUTPUT f 


(X) a THE SOLUTION, IT HAS NCOLSa ROwS AND NRHSS COLUMNS , 

CATNV) e THE LEAST SQUARE INVERS OF (A) AFTER THE REQUESTED 

ROWS and columns are deleted* (aInv) has ncOlSa rows 
and NROWSa columns, 

(RMSRES)a the SQUARE ROOTS OF THE AVERAGE SQUARED ERRORS FOR EACH 


r»f>r»r»r»or*r>r»r»r*or>r>nor»r>r>oooorir*rir»r»or»f> 


set OF kIGHT H4NO SIDES. 

DET I The determinant of thE least square matrix, 

lERR a ERROR CODE, 

0 means no error, 

1 MEANS READA»,TRUE, AND IFMT(WaS NOT 0,1, i, OR 3, 

? means there were fewer equations than unknowns 

(NROWSa,LT,NCOLSAI , 

X means nb less Than nrohSb 

OTHER ARGUMENTSf 


UTA) > storage space for performing calculations, ThE number 

OF ROWS IS NCOLSA AND THE NUMBER OF COLUMNS 
IS NCOLSA+NRHSS { I NVEHS# , F AL$E , ) OR NCOLSA* 
NRhS$*NR 0W3A (iNVEP8a,TRUF«) • 

(AR) m STORAGE SPACE, THE NUMBER OF ELEMENTS USED a THE 
number of columns of (ATA), 

(ID » STORAGE SPACE, THE NUMbER OF ELEMENTS USED a NrOLSA. 
OTHER notes* 


1, THERE IS A storage ADVANTAGE IN THIS PROGRaH SINCE 
only 1 ROW OF (A) NEEDS TO BE IN CORE aR ONE TIME. 
(INDCOL) will be COMPUTED AUTOMATICALLY |F NC0L5A«NC0LSB , 
5. (INDROm) will be computed AUTOMATICALLY IF NR0wSA«nR0WS8 . 
U, NB SHOULD be ,GE, NROWSb, 

IF READAa. FALSE, , TmEN NROWS ShOULD BE .r,E. NROwSfi, 


LOGICAL REA0a,INVER8 

DIMENSION a(NROW8iNCOL8B)#8(NB,NRhS8)# aTA(NCOLSa#NkTP) r 
I ARrNKTP),lLfNCOLSA)»lNftCOLf NCOLSA) #lNRROWeNROw$A)#X(NX,NRHSSl 
2A1 NV(Na1NV,NROWSa) >RMSRES(NRhSS) 

810 FORMATfftFlO.O) 

820 FORMATfSFlS.O) 

850 format (8f10«0) 


r> ri o o o r* 


C 

c ERROR CHECKS 

IF(REAOA)GO TO J 
NREcD « 1 
GO TO ? 

1 lP(lFMT,Lt,3,ANO,lFMT,GE,0) GO TO 2 
lERR « 1 
RETURN 

? IF{NR,GE,NH0WSB) go to 3 
IE«R « 3 
RETURN 

3 IF(NROwJ5a,GE,NCOLSa 3 GO TO S 
lE^R » 2 
return 

H IF fNCOtSA,NE,NCOLSe) GO TO T 
DO 6 I « 
h INOCOU(I) 9 I 

7 IE (NROwSa*NE,NRONSB) go to R 
DO 6 I * 1 ,nR0nSa 
g lNDROi«(I) 8 I 

INITIALIZATION 

clear ata array 

Q DO 10 J « I ,NCOLSA 
DO 10 K » 1 ,NkTP 
10 ATA(J,K) s 0,0 
TRSTOR e NRQNS 
IXX # IZERf) ♦ t 
JXX 8 IZeRO ♦ NRhSS 
IRBEGN s 1 
tPTROW 8 I 

IX s NCOLSA 4 NRhSS 4> 1 
INV » IX 

JFHTCK « IFMT ♦ 1 

read in one logical record or as many rows as you wish TO 
operate on at a time 

DO 200 K 8 1,NRECD 
lEf ,NOT,ReAOA) GO TO Ufi 


r» r» r» 


ICHICK « NR0NSB-fK-n*N«0W8 
UNO a MINo(NROWS,ICHECK) 

GO TO MOr20f2lf2?W IFMTcK 
C CONTROL SHOULD NgV^R TO HeRe 
STOP IS 

19 RE4D(NUNlT)(CA(L^N),N«i,NC0U5B),Ltl#ieND) 

GO TO ?6 

20 Rf AO(NUNTTf6lO)<(A(LiN)#Nil^NCnLSBl,Lal»IENoi 
GO TO ?b 

21 REAO(NUNlT^e2o)((ACLfN),Na|,NCOLSB>,L*i#UNOl 
60 TO ?b 

22 pEAOtNUNlTfSSO) ((A(L>N)^N»l»NCnLSB),L«W UND) 

C CHECK TO be sure RE WaNT TO REtAlN aT LEaST ONE ROW 

20 DO So J « IHBEGN.IRSTOP 

If (INOROW(lPTROW) , eO,J) go TO 9o 
SO CONTINUE 
GO TO IRO 

begin moving elements to AR array checking that This a row and 

COLUMN WE WISH TO RETAIN', 
an jROw » 1 

DO 170 J » IRBEGN, IRSTOP 
lFaPTROW,CT,N«OwSA) GO TO ITo 
IF(INDROW(IPTROW),ne,J) go to 170 

DO SO L 8 I^NCOlSB 
ie(h,gt ncolsa^ go to 50 
1F(1NDc5uH),Ne,L) go to So 

AR(M1 ■ A(lRow,L) 

M « H4 1 
50 CONTINUE 

I » INOROW(IPTROW) 

OO 60 L « IXX,JXX 
AR(H) s B(NL) 

60 M 8 M+l 

IPTROW « IPTROW + 1 
If C*N0 T,INVeR8) go To 75 
DO 70 L * jX,NKTP 
70 AR<L) ■ 0.0 


ARfJNV) 9 UO 
INV 8 INV ♦ I 

75 DP 90 I * l^NCOLSA 

R B SqRT(ATA(I,I)**2 ♦ arc!)**?) 

IP t»,EP#O,0) GO TO <>0 
C * ATA(I,1 )/r 
S * AR(n / « 

DO eo L « i,nktp 

T? • C • ATA(I,L) + S * ARfL) 

ARfl) « -S * ATAn»l) + C * ARfL) 

80 ATAfliL) • TJ? 

RO COWTlNUe 
170 IPOW * iROh ♦ 1 

iFfREAnA) GO TO too 
GO TO 205 

too JRBEGN 9 IRSTOP « \ 

IRSTOP m IR5T0P ♦ It^o 

200 continue 

205 II » t 

DO 220 1 * l.NCOLSA 
IF(ATA(I,n,LE.O,OOOOOOn GO TO 2l0 

11(1) * ir 
n 8 u 4 1 
GO TO 22o ‘ 

2lO IL(I) 9 0 
220 CONTINUE 

DO 365 vl » l,N«H$S 
NKTJ s NCOLSA ^ J 
DO 500 I « t, NCOLSA 
300 ARflt z 0,0 
II c NfOLSA 
DO 3«0 I 8 I, NCOLSA 
IF(lLCII),Lr,0) GO TO 3«0 
JI 8 IKH) 

IF(II^NCOLSA) 3i0f350»350 
3l 0 IF 9 I ! ♦ 1 

DO 320 K 8 IK, NCOLSA 


n r» o 


320 ARfll) s AR<!I) - ATA(JI»K) * AR(K) 

330 ARfll) » (ARfin ♦ ATAf JNNKTJ))/ATA(Jlf in 
3ao II • n • i 
350 CONTINUE 

DO 360 I ■ l,NCOL$A 
X(IfJ) • ARf!) 

360 CONTINUE 
565 CONTINUE 
OET « 1, 

DO 370 1 • 1,NC0lSA 
370 OfT « OET * ATAflf I) 

ir f .NOT, iNvgRSjGo TO aoo 

C 

c eino inverse 

DO 3R0 J « I .NCOUSA 
NKTj B NCOLSA 4 NRmSS ^ J 
DO 373 I • l.NCOtSA 
373 AR(I) • 0,0 
II ■ NcOlSA 
DO 360 I ■ 1 , NCOLSA 
lEflUin.LE.O.O) GO TO 380 

U1 » Ilfin 

IE flI*NC0L8A) 375, 379,585 

375 IK » 11 ♦ 1 

OO 377 K 8 IK, NCOLSA 
377 ARfll) 8 ARfll) . ATA(J1,K) * ARfK) 

579 AR(II) »(ARCI1) ♦ ATA(JI,NKTJ))/ATA(JIiII) 

380 II • 11 - I 
385 CONTINUE 

00 388 I « j, NCOLSA 
388 AinV( 1,J) 8 ARfj) 

590 CONTINUE 

aoo • 0 
Return 
end 

subroutine KRUNCHfNl, IND0EL^N2, INdEX,NLI) 

THIS subroutine COMPUTES THE LOGICAL UNION OF THE SET 

,,,,,0F integers in fINDEX) AND TmE SET OF INTErE^S NOT 
,,,,,1N (INddEL). 


oonor»r»nonr> 



THE INTEfiERS IN (INOOEUJ aNO (INDEX) A«£ ASSUMED TO hE 

.arranged In ascending ORDER N) aND Nj aRE TME 

OP EUEMENT5 IN (INODEU) a^C) (INDEX), RESREcT I VEUY . 
.....ON RETURN, The union is stored IN (INOeX) and Th£Re aRe 
.....N u elements. 

.... .example 

(JNDDEU) a 2r 6# 8» 1 0^ j 5 

.•••.AMO (INDEX) « Ii2»5,a,g,l0,ll,l2 , 

THEN ON RETURN ( I NDE X ) * 1 , 5# 1 1 , 1 ? AND NUi«, 

OlHENSION INooEUfNt) 
dimension INoEX(N?) 
n2Ti) 

N^HIcN^•i 

NUsn;> 

00 20 tal.Nl 
00 10 J»N2i,N2 
K».J 

ir(rNDEX(J),f,T,lND0EL(l)) GO TO J5 
IFfINDEXfJ),LT,lNDOfcL(l>) GO TO lO 
NUsNU*} 

1F(J,Ge,N2)ReTURN 
DO fi LaJ#N2M^ 
iNoex(L)»iNOf xd^n 
H CONTINUE 
00 TO i5 
10 CONTINUE 
l5 N2I»K 

IF(N?I,CT,N2)RfTURn 
20 CONTINUE 
Return 
End 

subroutine ST0P2(N,MESaGE#VaL) 

DIMENSION ME$*GE(2n) 

NRITE(NiI) MgSACEfVAl 
STOP 13 

1 FORHaT(///» ***** W20A^»f» *****»/» VAt«' ^ IPE IS,7) 

end 



